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13.  ABSTRACT  (Miximum  200worOI) 

The  use  of  new  advanced  solid  rocket  propellant  constituents,  such  as,  e.g., 
boron  as  HEDM  and  ammonium  nitrate  as  replacement  for  halogen-containing 
oxidants,  requires  knowledge  of  their  combustion  characteristics.  To  this  end 
experimental  rate  coefficient  k  measurements  on  individual  important  A1  and  B 
species  reactions  were  made  in  the  300-1800  K  temperature  range  using  a  high- 
temperature  fast-flow  reactor.  Temperature  dependent  rate  coefficient  expressions^ 
km  =  AT(super  n)exp(-E/RT)  are  found  for  the  reactions  between  AlO  and  CH(sub  4), 
0(sub  2);  A1  and  N(sub  2)0;  A1C1  and  N(sub  2)0,  CO(sub  2),  SO(sub  2);  BC1  and  N(sub 
2)0,  SO(sub  2);  BO  and  CO(sub  2).  To  allow  semi-empirical  predictions  for  further 
metal  oxidation  reactions,  two  semi-empirical  theories  have  been  developed  and 
utilized  These  are  Semi-Empirical  Configuration  Interaction  theory  for  prediction  of  E 
from  one  species  to  another  and  Metals-Transition  State  Theory  for  prediction  of  k(T) 
from  k  at  one  temperature.  Dual  military-civilian  applications  of  this  work  are 

discussed. 
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I.  RESEARCH  OBJECTIVES 


Advances  in  solid-fuel  propellant  technology  for  the  late  nineteen  nineties 
depend  critically  on  an  understanding  and  quantitative  description  of  a  number  of 
combustion  processes.  These  include  the  combustion  of  metallic,  particularly  boron, 
HEDMs  and  environmentally  clean  A1  and  B  propellants,  containing,  e.g.  ammonium 
nitrate  instead  of  perchlorate.  In  this  context  the  present  work  had  several 
objectives.  Unique  high-temperature,  fast-flow  reactors  HTFFR,  which  allow 
measurements  of  the  kinetics  of  isolated  reactions  in  the  300-1900  K  temperature 
range,  were  employed  to  provide  a  kinetic  data  base  for  metallic  species  reactions  in 
propulsion  system  models.  The  results  obtained  enhance  the  understanding  of  the 
chemistry  involved,  which  in  turn  allowed  making  predictions  for  further  reactions. 
This  enables  faster  model  development  than  experimental  work  by  itself  can  provide. 
To  this  end  establishing  correlations  in  series  of  similar  reactions  and  understanding 
the  fundamental  basis  for  these  was  emphasized.  Where  appropriate  dual  military- 
civilian  use  aspects  of  the  work  have  been  emphasized  (see  especially  Section  III  A). 

n.  RESULTS 
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(1993). 
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HI.  SIGNIFICANCE  OF  THE  RESULTS 

Here  we  discuss  the  significance  of  the  papers  given  in  the  preceding  section  in 
the  chronological  sequence  given  there. 

A.  General  and  Military-Civilian  Dual-Use  (Publications  1  and  2) 

The  general  significance  of  Publication  1,  should  be  clear  from  its  Preface  given. 
It  has  become  very  widely  quoted. 

A  few  dual  military-civilian  aspects  are  illustrated  here.  Publication  2  is  one 
example  of  this;  another  is  the  use  of  our  unique  family  of  reactors  developed  for 
wide  temperature  range  studies  on  isolated  reactions.  The  development  of  these 
reactors  was  prompted  by  the  needs  of  the  military  for  kinetic  data  on  metal 
combustion.  With  support  of  NSF  we  are  now  also  applying  the  reactors  to  waste 
incineration  problems.  In  that  area  we  showed  the  presence  of  several  paths  for 
chromium  oxidation.1  As  chromium  is  in  its  highest  oxidation  state  highly  toxic,  but 
relatively  harmless  in  lower  oxidation  states,  knowledge  of  the  various  reaction 
pathways  is  needed  to  hinder  the  complete  oxidation.  An  understanding  of  high- 
temperature  reactions  of  metals  is  important  to  many  more  applications.  These 
include  decreasing  slag  formation  in  coal  combustion,  preventing  hot  corrosion  in  jet 
engine  turbines,  alleviation  of  soot  formation  in  all  sorts  of  combustion  devices,  as 
well  as  new  chemical  synthesis  processes  for  refractory  materials. 

The  direct  conversion  of  methane,  i.e.  natural  gas,  to  liquid  fuels  and  chemicals 
of  commercial  interest  is  a  long- sought  goal.  A  major  approach  being  followed  is  to 
pass  methane  over  solid-state  catalysts  consisting  of  mixtures  of  metal  oxides  with 
dopants.  The  first  step  in  this  process  is  to  abstract  a  hydrogen  atom  from  methane 
CH4  to  produce  a  highly  reactive  methyl  radical  CH3  to  form  C-C  bonds.  To 
understand  these  surface  reactions  better,  gas-phase  studies  of  the  reactions  between 
methane  and  metal  oxides  are  desired.  Theoretical  studies  have  suggested  that  such 


2 


abstraction  reactions  should  occur  at  elevated  temperatures.  Because  of  our 
experience  with  AlO  from  studies  for  rocket  propellants,  we  decided  to  use  that  oxide 
to  experimentally  investigate  this.  Publication  2  provides  the  first  experimental 
confirmation,  but  also  shows  that  at  least  in  this  instance  the  theoretically  predicted 
mechanism  is  not  dominant.  In  general,  such  gas-phase  studies  provide  essential 
input  for  gas-solid  catalytic  reactions. 

B.  k(T)  Predictions  for  Model  Development  (Publications  3-5) 

Combustion  models,  whether  for  the  plume  or  chamber,  require  reliable  kinetic 
input  data  at  realistic  temperatures.  The  data  needed  are  rate  coefficients  k(T)  and 
mechanisms.  Measurements  with  our  HTFFR  (high-temperature  fast-flow  reactor) 
and  MHTP  (metals  high-temperature  photochemistry)  techniques  are  providing  such 
information  for  metallic  species.  However,  experiments  alone  are  a  rather  slow  way 
to  obtain  information  on  all  the  reactions  needed.  It  is  therefore  desirable  to  develop 
predictive  techniques.  Fully  theoretical  a  priori  techniques  have  not  yet  reached  the 
state  of  making  meaningful  predictions.  We  have  therefore  chosen  to  develop  semi- 
empirical  techniques  for  metallic  species.  These  predictions  allow  either  extrapolation 
of  our  measured  data  from  one  species  to  another  (SECI  theory),  or  for  one  reaction, 
from  one  temperature  to  a  wide  range  of  temperatures  (MTST  theory).  Another  way 
of  expressing  the  scope  of  these  theories  is  that  in  rate  coefficient  expressions 

k(T)  =  AT11  exp(-E/RT)  (1) 

SECI  predicts  the  exponential,  MTST  the  pre-exponential  part.  The  theories  are  semi- 
empirical  in  that  for  SECI  the  E  factor  (activation  barrier)  of  a  similar  reaction  needs 
to  be  known  and  that  for  MTST  an  experimental  rate  coefficient  measurement  at  one 
temperature  is  required. 

The  left  hand  side  of  Figure  1  illustrates  SECI  the  right  hand  side  MTST,  each  of 
which  will  now  be  discussed. 

1 .  Semi-empirical  Configuration  Interaction  (SECI)  Theory 

The  E  factor  in  eq.  (1)  is  the  difference  in  the  potential  energy  of  the  activated 
complex  (transition  state),  formed  from  the  reactants,  and  the  sum  of  the  potential 
energies  of  the  reactants.  It  represents  the  highest  potential  energy  on  the  path  from 
reactants  to  products.  It  can  be  calculated,  using  resonance  theory,  from  the 
properties  of  three  hypothetical  resonating  activated  complexes,  i.e.  that  resulting 
from  purely  ground  state  interaction  for  which  the  barrier  is  q,  the  interaction 
involving  the  lowest  excited  state  of  the  atomic  reaction  partner,  and  ionic  interaction. 
The  wave  function  of  the  actual  complex  resulting  from  these  structures  is  given  by 

*¥  =  Cl'Fi  +  C2vF2  +  C3W3  (2) 

where  the  c  factors  give  the  contributions  of  each  structure.  The  actual  'F  will  be  that 

for  which  E  =  J  vF*H'Fdx  is  at  a  minimum.  The  condition  can  be  expressed  as  a  set  of 
simultaneous  linear  equations,  ultimately  resulting  in 
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(3) 


aq  +  bq2  +  cqE  +  dq2E  +  fE  +  gE2  +  hE^  —  0 

where  a,  b,  c,  d,  e,  f,  g,  h  are  constants  derived  from  tabulated  physical  properties  of 
the  individual  reaction  partners. 

E  and  q  can  not  simultaneously  be  calculated  from  first  principles.  For  a  series 
of  similar  reactions,  e.g.  the  boron  group  atoms  with  N2O  (see  the  figure)  we  therefore 
obtain  E  for  one  atomic  partner.  This  is  done  by  fixing  n  for  the  series,  which  allows 
finding  E  from  eq.  (1)  used  as  the  fitting  expression  for  the  experimental  k(T) 
measurements.  Eq.  (3)  then  yields  q,  which  according  to  the  Hirschfelder  theory  is 
the  same  for  all  members  of  the  series.  E  for  the  other  members  of  the  series  can 
now  be  calculated  from  eq.  3  using  q  and  the  constants  appropriate  for  each  of  the 
atoms. 

The  SECI  figure  shows  E  plotted  against  I.P.  (the  ionization  potential  related  to 
the  3rd  structure)  plus  P.E.  (the  promotion,  or  excitation,  energy  representative  of  the 
2nd  structure).  It  shows  good  correlation  with  IP  +  PE  and  excellent  agreement 
between  theory  and  experiment  obtained  for  three  groups  of  atom  reactions  with  N2O. 
The  left  hand  line  is  for  the  one-electron  alkali  metal  group  and  is  of  interest  for 
example  to  hot  turbine  combustion  corrosion  problems.  The  middle  line  is  for  two 
electron  atoms,  and  contains  metals  of  interest,  e.g.  to  toxic  emission  problems.  The 
right  hand  line,  for  two  s  one  p  electron  atoms,  has  most  recently  been  obtained  and 
includes  metals  of  most  importance  to  rocket  propulsion. 

2.  Metals  Transition  State  Theory  (MTST) 

As  shown  in  Fig.  1,  by  taking  one  experimental  point  (the  solid  symbol)  at  mid- 
temperature  range  for  the  reaction  of  interest,  we  can  predict  the  shape  of  the 
Arrhenius  plot  in  good  agreement  with  measurements,  thus  reducing  the  work 
required  for  data  taking. 

In  general  transition  state  theory  (TST)  k(T)  is  given  by: 

k(T)  =  q;~  exp(-AE*/RT)  (4) 

'  h  qABqCD 

where  ke  and  h  are,  respectively,  the  Boltzmann  and  Planck  constants,  and  qAB  and 
qCD  are  the  partition  functions  of  the  reactants,  q*  and  AE*  are,  respectively,  the 
partition  function  of  the  transition  state  and  its  energy  above  the  reactant  ground 
states,  both  of  which  are  unknown.  By  using  simple  ab  initio  methods  for  estimating 
q*,  and  equating  eq.  (4)  with  an  experimental  rate  coefficient  at  one  temperature,  AE* 
can  be  found. 

The  basic  TST  approach  was  developed  primarily  for  reactions  of  more 
commonly  studied  reactants,  such  as  hydrocarbons.  There,  the  transition  state 
structure  and  hence  q*  can  be  found  readily  by  comparison  to  known  stable 
compounds,  which  information  is  generally  lacking  for  refractory  species.  The 
advance  is  that  we  have  made  the  method  useful  for  such  species  as  well. 
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C.  AIO  and  BO  Oxidation  Reactions  (Publication  6) 


Oxidation  reactions  of  A1  and  B  species  are  essential  input  to  various  rocket 
combustion  models.  In  further  experimental  work,  we  have  thoroughly  studied  the 
AIO  oxidation  by  O2  (Publication  6)  and  obtained  a  significant  preliminary  result  for 
BO  oxidation  (not  published). 

The  AIO  +  O2  reaction  system  was  originally  thought  to  proceed  as  a  simple 
abstraction  reaction  leading  to  AIO2  +  O.2  However,  in  recent  years  we  have  found 
several  examples  of  reactions,  considered  in  the  literature  as  simple  abstraction  to 
proceed  by  addition  or  insertion  at  lower  temperatures,  with  abstraction  becoming 
more  significant,  but  not  necessarily  dominant,  at  higher  temperatures.  Examples  of 
this  are  publication  2  and  Ref.  1.  These  observations  and  the  vast  improvements  in 
HTFFR  technique  since  then,2  prompted  this  new  investigation.  For  temperatures 
below  about  1000  K  the  AIO  consumption  now  is  found  to  be  pressure-dependent, 
with  a  negative  temperature  dependence;  log  k  (305-1010  K)  =  -25.36-1.69  log  (T/K) 
cm6molecule-2s-1.  At  higher  temperatures  the  system  becomes  pressure- 
independent  and  the  AIO  consumption  rate  coefficients  increase;  k(  1200- 1690  K)  = 
7.7xl0*10exp  (-10008  K/T)  cm3molecule*1s-1. 

The  following  mechanism  describes  these  findings: 


AlO  +  O- 
I  1 
I 


aio2+o 


AIO3 


(5) 


Here,  the  dotted  lines  indicate  the  paths  unlikely  to  contribute  significantly.  Thus,  the 
reaction  is  thought  to  proceed  through  an  intermediate  AIO3  complex,  which  can  be 
collision-stabilized  or  dissociate  to  AIO2  +  O.  The  latter  process  is  96  +  40  kJ  mol*1 
endothermic. 


It  is  interesting  to  observe  that  in  a  room  temperature  study  of  the  AIO  +  CO2 
reaction,  competition  between  an  addition  and  abstraction  reaction  was  observed.3 
This  suggests  that  in  that  system  abstraction  will  become  dominant  at  lower 
temperatures  than  for  AIO  +  O2.  The  O2  reaction  should  be  primarily  of  interest  for 
plume  models,  while  the  CO2  reaction  is  of  both  plume  and  chamber  interest. 

The  AIO  +  O2  behavior  contrasts  with  that  of  the  iso-electronic  exothermic 
BO  +  O2  reaction  system,  for  which  no  pressure  dependence  was  found  over  the  300- 
1000  K  range,  but  again  a  negative  temperature  dependence  was  observed;4 
k(T)  =7.0xl0*12exp  (259  K/T)  cm3molecule*1s*1.  This  suggests  a 

BO  +  O2  ^  BO3  -4BO2  +  O  (6) 

mechanism.  We  have  made  a  series  of  BO  oxidation  measurements  in  the  HTFFR 
apparatus.  In  this  environment  accurate  rate  coefficients  could  not  be  obtained, 
apparently  due  to  BO  wall-dimerization.  We  are  about  to  continue  this  work  in  one  of 
our  MHTP  (metals  high-temperature  photochemistry)  reactors,  where  such  problems 
should  not  arise.6  Nonetheless,  the  HTFFR  work  has  already  yielded  some  interesting 
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information.  Our  results  on  BO  +  O2  between  400  and  1000  K  are  in  essential 
agreement  with  the  earlier  measurements  from  an  MHTP-like  reactor.4  For  the 
hypothesized  reaction  BO  +  CO2  -»  BO2  +  CO  a  slow  pressure-dependent  process  was 
observed  at  room  temperature,  but  no  measurable  BO  consumption  occurred  from  780 
to  1200  K.  This  indicates  rate  coefficients  less  than  about  lxl0‘15  cm3molecule-1s*1. 
This  is  four  orders  of  magnitude  lower  than  found  for  the  BO  +  O2  reaction.  It  further 
is  two  to  three  orders  of  magnitude  lower  than  the  guestimated  values  in  the  current 
models  for  boron-hydrocarbon  slurry  (air-breathing)  combustion,  where  it  is  indicated 
to  be  a  highly  sensitive  reaction.  Those  models  should  now  be  reconsidered. 
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a  variety  of  basic,  technical arable  number  oLscientists  and 
currently  developed  world  metal  reactions.  Yet.  there  has  been 

engineers  occupied  with  ga  P  rentiy  n0  book  or  major  meeting  has 
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previously  been  dedicated  ,  .  .  uc  during  a  recent  sabbatical 

With  colleagues  our  knowledge.  The  present 

clearly  brought  out  the  n  Qnse  t0  thli  it  brings  together  overviews 

volume  has  come  about  s  ’  wuh  fundamental  research,  both  of  an 

from  most  of  the  gr  P*  .  nature  ^till  much  of  this  work  has  been 
experimental  and  theoret  ca  m  Th  individual  authors  show  these 

stimulated  by  very  practical l  prob  m,  ^ne  &  few  of  the  final  chapters 
connections.  To  further  emp 'has  what  research  needs  are.  or 

primarily  address  their  solution.  The 

how  fundamental  informa  connection  between  the  various 

'"‘her  eff»rts  and  a,,pl,cauons- 
if  one  compares  the  present  field  to 

studies  such  »  .  al  of  the  theory  chapters  in  this  volume 

effort  there  has  been.  In  ta  • -  ,  ,  ror  u  stimulated  the  research, 

came  about  primarily  becaus  between  the  kinetic  and  dynamic 

An0,h,'r  £&%£*»£  dVn«  Wo™  a  proper  synthesis  of  the  results 
aspects,  much  work  neea  he  chapter5  indicate  current 

obtained  can  be  achieved.  -  important  contnbutions  can  be 

further  r'esearch  and  interdisciplinary  collaboi  ation. 

The  birth  of  this  volume  has  proceeded  £  emed  and  discussed 

submitted  preliminary  'j.tiaPt^rs  j  chemistrv  Division  of  the  American 

at  a  symposium  held  oetore  the  Following  the  meeting. 

Chemical  society,  in  New  Wk  irV  “u  were  then  refereed  and  usually 

»-  .ho  benefit  of  much 

interaction  on  an  international  scale 
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A  recent  ab  inito  study  suggested  that  gaseous  metal  oxides  can  directly  abstract  H  atoms  from  hydrocarbons, 
but  with  considerable  barriers.  To  try  to  confirm  the  occurrence  of  such  reactions,  the  title  reaction  has  ten 
studied  in  a  high-temperature  fast-flow  reactor.  The  data  are  well  described  by  the  fitting  expression  *(590- 
H80  IQ  ■  9  6  X  10-”(77K)7  %  exp(2468  K/D  cm3  molecule'1  s*1,  with  2<r  precision  limits  varying  with 
temperature  from  ±4%  to  ±12%.  The  corresponding  2<r  accuracy  limits  are  about  ±25%.  C°m^n»n  ofth« 
result  to  further  ab  initio  and  conventional  transition-state  theory  calculations  suggests  that  direct  abstraction 
indeed  can  occur  at  the  hi^ier  temperatures  but  that  one  or  more  other  channels,  possibly  involving  AlO 
insertion  into  a  C-H  bond,  dominate  in  the  initial  attack  step.  The  potential  implications  of  this  work  for 
catalytic  conversion  of  methane  to  higher  hydrocarbons  are  considered. 


Introduction 

In  a  recent  ab  initio  study  Borve  and  Pettersson1  predicted 
that  some  gaseous  metal  oxides,  LiO,  MgO,  and  AlO,  can  directly 
abstract  an  H  atom  from  methane  to  produce  methyl  radicals. 
Such  exothermic  reactions  have  not  been  observed  in  the  gas 
phase,  and  their  studies  suggested  these  to  have  considerable 
barriers,  25—67  kJ  mol'1 .  However,  solid  metal  oxides  have  been 
found  to  be  effective  catalysts  for  H  abstraction  from  methane, 
leading  to  C2H4  and  C2H«  in  significant  yields,2  a  process  of 
considerable  commercial  promise.  The  basis  for  this  catalytic 
reaction  has  been  strongly  suggested  to  be  the  surface  reaction 
of  the  O"  ion  to  remove  a  hydrogen  atom,  a  process  which  appears 
to  proceed  to  a  significant  degree  only  above  about  1000  K.3 
Theoretical  studies  of  cluster  models  of  the  M0O3  oxide  structure 
have  supported  this  hypothesis.4  If  the  occurrence  of  the  gas- 
phase  processes  could  be  experimentally  confirmed,  it  would  not 
only  establish  a  new  class  of  reactions  but  also  could  lead  to  a 
better  understanding  of  the  heterogeneous  catalytic  processes. 
We  therefore  decided  to  investigate  one  of  these  suggested 
reactions: 

AlO  +  CH4  —  AlOH  +  CH3  (1) 

for  which  A//0  =  -27  kJ  mol-'.s-6  While  this  reaction  has  a 
higher  predicted  barrier  than  those  of  the  other  oxides,1  we  have 
had  extensive  experience  with  the  study  of  the  gas-phase  reactions 
of  AlO  at  temperatures  up  to  1600  K  in  HTFFRs  (high- 
temperature  fast-flow  reactors)7  and  therefore  selected  this  oxide. 

Technique 

The  HTFFR  facility  used  in  this  work  and  the  data  handling 
procedures  have  been  previously  described/  Briefly,  a  vertical 
mullite  (McDanel  MV-30)  reaction  tube  (60  cm  long,  2.2  cm 
i.d.)  is  heated  to  the  desired  temperature  by  columns  of  resistively 
heated  SiC  rods  inside  an  insulated,  water-cooled  vacuum  housing. 
A1  vapor  is  produced  by  resistively  heating  an  Al-wetted  W  coil 
in  a  flow  of  Ar  carrier  gas.  A  trace  of  oxidizer,  0.5%  N20  or  5% 
02  in  Ar,  at  X  10~3%  of  the  mam  Ar  flow,  is  passed  through 
a  Pt  side  tube  to  a  location  just  downstream  from  the  vapor  source 
to  rapidly  convert  A1  to  AlO.8-10  Further  downstream  CH4/ Ar 
mixtures  at  flow  rates  of  4—8%  of  the  main  Ar  flow  rate  are 
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introduced  through  a  movable  inlet.  Rate  coefficient  measure¬ 
ments  are  made  under  pseudo-first-order  conditions,  [AlO]  « 
[CH4],  in  the  stationary  inlet  mode,11  at  reaction  zone  lengths 
of  10  or  20  cm. 

Relative  AlO  concentrations  were  monitored  by  laser-induced 
fluorescence  (LIF)  using  a  pulsed  Lambda  Physik  EMG  101 
excimer/FL  2002  dye  laser  in  combination  with  a  KDP  doubling 
crystal.  Two  transitions  B22-X22  and  C22-X22  were  used.  In 
the  former,  AlO  is  pumped  on  the  464.8-nm  (1,0)  band  and  the 
fluorescence  is  observed  through  a  482-nm  (20-nm  fwhm) 
interference  filter,  i.e.,  mainly  at  the  486.6-nm  (1,1)  band. !  2  The 
B2S-X22  system  could  only  be  used  up  to  1200  K  because  of 
interference  by  the  background  radiation  from  the  reactor  walls. 
For  measurements  above  1200  K  and  some  checks  at  lower 
temperatures,  the  C2Z-X22  (0,0)  transition  at  302.2  nm12  was 
used  in  combination  with  a  30l-nm  (1 1-nm  fwhm)  interference 
filter.  The  fluorescence  was  detected  by  an  EMI  9813QA 
photomultiplier  tube,  connected  to  a  Data  Precision  Analogic 
6000/620  100-MHz  transient  digitizer. 

The  gases  used  were  99.998%  Ar  from  the  liquid  Ar  (Linde), 
0.5%  N:0  (99.99%)  in  Ar  (99.999%)  from  Matheson,  5%  02 
(99.99%)  in  Ar  (99.995%)  from  Scott,  and  CH4  (99.99%)  from 
Matheson. 


Results 

Plots  of  ln[A10)reuu«c  versus  [CH.]  yield  straight  lines  with 
slopes  -kt,  where  t  is  reaction  time.  For  each  individual 
measurement,  k  and  ak  were  determined  by  a  weighted  linear 
regression.1113  The  k  values  and  the  experimental  conditions 
under  which  they  were  obtained  are  summarized  in  Table  I.  The 
measurements  may  be  seen  to  be  independent  of  the  following: 
pressure  varied  from  13.1  to  63.3  mbar,  corresponding  to  total 
concentrations  [M]  from  8.0  X  10“  to  6.4  x  1011  molecules 
cm-3;  average  gas  velocity  v  varied  from  12  to  98  m  s  1 ,  reaction 
zone  length  selected  at  10  or  20  cm;  initial  fluorescence  intensity 
F  (a  measure  of  initial  [AlO])  varied  from  9  to  65  in  arbitrary 

units.  The  temperature  range  covered  is  590-1 380  K.  The  lower 
limit  was  determined  by  the  heating  effect  of  the  A1  vaporizer, 
and  the  upper  limit  by  the  dissociation  of  CH..  An  Arrhenius 
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TABLE  Is 


Summary  of  R»te  Coefficient  Measurements  on  the  AlO  +  CH,  Reaction' 


reaction  zone 
length  (cm) 


p  (mbar)  [M]  ( 10'7  cm  •')  [CH4]ma*  ( 10' 5  cm  ' ) 


V(ms')  T(  K) 


k  ±  <n  (cm3  molecule  1  s  l) 

4.42  ±  C.33  (-14)'' 
4.93  ±  0.41  (-14)' 

3.19  ±0.31  (-14)' 

5.22  ±0.63  (-14)' 

6.78  ±0.93  (-14)' 

1.67  ±0.21  (-13)' 

1.83  ±0.21  (-13)' 

1.76  ±0.23  (-13)4 
1.72  ±0.21  (-13)' 

1.75  ±0.17  (-13)' 

1.09  ±0.12  (-13)' 

1.58  ±0.19  (-13)' 

1.81  ±0.22  (-13)' 
3.30  ±0.53  (-14)' 

8.15  ±0.90  (-14)' 
8.53  ±0.79  (-14)' 

'  1.11  ±0.08  (-13)' 

8.82  ±0.97  (-14)' 
8.62  ±0.81  (-14)' 
1.43  ±0.16  (-13)' 

9.28  ±  1.07  (-14)' 

1.28  ±0.14  (-13)' 

1.29  ±0.15  (-13)' 
7.38  ±0.53  (-14)' 

7.22  ±0.55  (-14)' 
7.91  ±0.69  (-14)' 

8.18  ±0.66  (-14)' 

7.11  ±0.63  (-14)' 

8.19  ±0.85  (-14)' 
1.09  ±0.09  (-13)' 
1.09  ±0.09  (-13)' 

2.10  ±0.21  (-14)' 

1.16  ±0.15  (-14)' 
1.34  ±0.14  (-14)' 
1.31  ±0.12  (-14)' 
1.65  ±0.22  (-14)' 
2.67  ±  0.32  (-14)' 

3.12  ±0.39  (-14)' 
7.05  ±0.83  (-14)' 
7.84  ±0.97  (-14)' 

1.10  ±0.12  (-13)' 
1.07  ±0.14  (-13)' 
1.44  ±0.18  (-14)' 
1.62  ±0.20  (-14)' 
1.99  ±0.14  (-14)' 
5.60  ±0.64  (-15)' 
9.95  ±  1.27  (-15)' 
1.02  ±0.13  (-14)' 

2.11  ±0.23  (-14)' 
7.78  ±0.91  (-15)' 

8.11  ±0.65  (-15)' 

6.71  ±0.77  (-15)' 
6.23  ±0.78  (-15)' 

1.72  ±0.20  (-14)' 
2.40  ±0.26  (-14)' 

2.43  ±  0.31  (-14)' 
4.94  ±0.54  (-14)' 
3.02  ±0.31  (-14)' 
4.64  ±0.53  (-14)' 

3.87  ±0.47  (-14)' 
5.81  ±  0.75  (-14)' 
4.57  ±0.32  (-13)' 
3.49  ±0.34  (-13)' 
1.46  ±0.11  (-14)/ 
1.36  ±0.10  (-14/ 
1.71  ±0.11  (-14/ 
1.68  ±0.14  (-14/ 
1.45  ±0.11  (-14/ 

3.44  ±0.37  (-13)' 
5.07  ±  0.44  (-13)' 

5.87  ±0.50  (-13)' 
4.52  ±0.35  (-13)' 
5.21  ±0.44  (-13)' 


20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

10 

10 

10 

10 
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20 

10 

10 

10 

10 

10 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

10 

20 

20 

20 

20 

20 

20 

20 

20 

10 

10 

20 

20 


26.7 

2.1 

26.5 

2.1 

43.3 

3.3 

20.0 

1.5 

20.0 

1.5 

15.5  ' 

1.0 

15.5 

10 

13.1 

0.8 

17.2 

1.1 

17.2 

1.1 

28.9 

1.9 

19.2 

1.2 

16.7 

1.1 

20.1 

1.6 

36.7 

2.6 

30.4 

2.1 

32.8 

2.2 

32.8 

2.3 

36.3 

2.5 

19.6 

1.4 

19.6 

1.4 

19.6 

1.3 

19.6 

1.3 

31.7 

2.3 

33.6 

2.5 

34.3 

2.6 

34.3 

2.6 

34.3 

2.6 

34.3 

2.6 

34.4 

2.5 

34.4 

2.5 

24.1 

2.1 

24.9 

2.5 

24.9 

2.5 

25.1 

2.5 

25.7 

2.5 

36.0 

3.3 

35.7 

3.3 

16.0 

1.2 

16.0 

1.2 

16.3 

1.3 

16.3 

1.3 

21.2 

2.4 

39.5 

4.3 

39.5 

4.3 

34.1 

4.1 

34.0 

4.0 

34.0 

3.9 

53.1 

5.3 

24.5 

2.9 

24.5 

2.9 

14.3 

1.7 

14.3 

1.7 

63.3 

6.4 

15.3 

1.4 

15.3 

1.4 

21.9 

1.8 

56.0 

4.8 

13.7 

1.2 

13.7 

1.2 

14.1 

1.2 

44.5 

2.3 

44.7 

2.5 

41.7 

4.5 

41.6 

4.6 

53.9 

5.3 

54.3 

5.6 

54.3 

5.7 

17.5 

1.0 

24.1 

1.4 

24.1 

1.3 

24.0 

1.3 

24.0 

1.3 

9.4 

9.4 
10.3 

4.8 

4.8 
5.0 
5.0 

4.4 

3.8 
3.8 

4.7 

5.8 

4.9 

9.9 

4.4 
5.8 

5.1 

8.5 

9.4 

6.4 

7.7 
5.0 
5.0 
7.0 

7.5 

9.6 
9.5 

7.4 

7.4 

5.7 

5.7 

13.7 

18.4 

18.1 
18.6 

18.8 

15.5 
15.2 

5.6 

4.5 

6.9 
8.1 

14.0 

17.1 

17.2 
21.1 

17.1 

16.9 

6.2 

20.6 
20.0 
11.0 

11.1 

9.9 

8.7 

8.8 

5.9 
4.7 
4.7 

4.7 
7.0 
1.1 
1.2 

7.4 

7.5 

6.8 

8.3 

8.5 

1.5 
2.0 
2.0 

1.3 
1.3 


17 

17 

24 

33 
24 
22 
20 
30 

47 
32 
38 
41 
64 

53 
30 
41 

38 
35 
32 
64 
63 

54 
53 

48 

34 

37 

39 
29 

24 

25 

23 

26 
28 

24 
28 
28 

52 

38 

53 

35 

39 
65 
34 

9 

9 

42 

30 

25 
9 

15 
27 

33 
20 
11 

34 

40 
25 

23 
20 

16 

24 

19 

20 
32 
45 
62 
65 
58 
40 
37 

37 

38 

39 


32 

33 
29 
62 
62 
98 
98 
90 
68 
68 
48 
56 
67 
46 

29 

33 
3! 

30 
28 
51 

51 

52 
52 

37 
35 

34 
34 
44 
44 
46 
46 

38 
29 
29 
28 
28 

34 

35 
58 
58 
57 
57 
28 
15 
15 
19 
19 
19 
21 

19 

20 
24 
23 
20 
37 
37 
33 
20 

55 

56 
56 
30 
28 

13 
12 

14 
13 
13 
60 

46 

47 

48 
48 


913 

912 

928 

932 

935 

1083 

1083 

1135 

1098 

1102 

1103 

1139 

1144 

882 

1001 

1051 

1053 

1031 

1034 

1034 

1037 

1061 

1063 

985 

982 

961 
963 

962 
961 
1000 
1003 

827 

721 

734 

717 

734 

781 

795 

945 

953 

938 

942 

651 

660 

663 

602 

618 

628 

727 

603 

621 

598 

591 

718 

780 

770 

893 

845 

846 
853 
866 

1374 
1380 
668 
650 
736 
707 
686 
1232 
1283 
1296 
•  1317 
1324 


-  The  measurements  are  reported  ,n  the  sequence  in  wh.ch  they  were  obtamed.  »  In  arbitrary  /  Dau 

obtained  by  using  the  C:2-X:I  transition  and  an  alumina  ring  inlet.  '  Data  obtained  by  using  the  C-2  X  2  transition  q  g 

obtained  by  using  the  B :2-X:2  transition  and  a  quartz  ring  inlet. 
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T,  K 


Figure  1. 

reaction: 


Arrhenius  plot  of  the  rate  coefficient  data  of  the  AlO  +  CH. 
O,  5%  O:  used  as  AlO  precursor;  0, 0.5%  NjO  used  as  AlO 


precursor. 


plot  of  these  data.  Figure  1,  shows  strong  curvature.  Several 
different  expressions  were  tried  to  fit  the  In  *  .ersus  T~'  data. 
Using  Marquardt’s  method,14  a  regression  fit  to  a  k(T)  =  A(T/ 
K)"  cxp(-B/T)  expression  yields 


*(590-1380  K)  =  9.6  x  10'3,(7yK)7%  X 

exp(2468  K/T)  cm3  molecule  1  s'1  (2) 

Equation  2  must  be  seen  as  a  fitting  expression  only,  as  the  numbers 
are  unrealistic  for  a  single  channel  reaction.  This  k(T)  expression 
most  likely  indicates  competition  between  several  reaction 
channels.  We  next  attempted  a  similar  regression  for  a  double¬ 
exponential  expression  k(T)  =  A  exp (-B/7)  +  C  exp (-ZJ/T), 
which  yielded 


*(590-1380  K)  =  2.0  X  10*10  exp(-8299  K/T)  + 

4.2  X  lO*13  exp(-2402  K/T)  cm3  molecule*1  s'1  (3) 

This  Fitting  expression  also  does  not  seem  to  be  representative  of 
individual  actual  reaction  paths,  because  of  the  high  preexpo¬ 
nential  of  the  first  term.  This  term  appears  too  large  for  a 
metathesis  reaction  between  a  diatomic  radical  and  a  molecule. 
While  molecular  insertion  reactions,  e.g.,  of  singlet  CH:  and  S1H2, 
are  known  with  A  factors  of  this  magnitude,  these  reactions 
proceed  with  an  activation  energy  close  to  zero.16-18 

Calculating  the  uncertainties  of  these  fitting  expressions  by 
combining19  variances  and  covariances,  similarly  as  in  our  earlier 
studies,7  yields  precision  limits  of  ±4%  to  ±12%  for  eq  2  and  of 
±5%  to  ±12%  for  eq  3,  depending  upon  temperature.  Both 
expressions  lead  to  essentially  equal  accuracy  limits  of  ±23%  to 
±26%,  when  combined  with  a  ±10%  uncertainty  in  the  reactor 
flow  profile11-20  and  ±20%  for  potential  systematic  errors.  The 
simpler  expression  (2)  is  therefore  preferred  for  the  present 
temperature  range. 


Discussion 

There  has  been  one  previous  attempt  to  observe  reaction  1 .  In 
that  room  temperature  study  no  reaction  could  be  detected,  and 
an  upper  limit  for  the  rate  coefficient  of  5  x  10" 14  cm3  molecule 
s-i  was  derived.8  Equations  2  and  3  suggest  much  lower  values; 
hence,  the  results  do  not  contradict  one  another. 

The  activation  energies,  i.e.,  the  local  slopes  of  the  Arrhenius 
plot,  Figure  1,  reach  the  67  kJ  mol-1  calculated  by  Borve  and 
Pettersson1  only  at  the  highest  temperatures  investigated.  For 
example,  from  eq  2  these  energies  are  46,  59, 66,  and  71  kJ  mol"1 
at  1000, 1200, 1300,  and  1380  K,  respectively.  Also,  the  fitting 
expressions  suggest  that  the  reaction  follows  several  paths. 
Further  mechanistic  considerations  thus  appear  in  order.  To  this 
end  we  first  assess  the  likely  importance  of  H-atom  abstraction. 


Berve  and  Pettersson  carried  out  a  partial  C  ASSCF  optimization 
of  an  assumed  C3t>  transition  state  (TS)  but  did  not  obtain  enough 
information  to  predict  rate  coefficients.  In  order  to  estimate 
these,  we  have  analyzed  abstraction  by  conventional  transition- 
state  theory  (CTST),  which  requires  vibrational  and  rotational 
data  to  estimate  the  partition  function  of  the  TS.  Accordingly, 
we  reoptimized  the  complete  TS  structure,  at  the  MP2/3- 
21G(*)  level,21  and  obtained  the  vibrational  frequencies.  The 
geometry  is  shown  in  Figure  2a  and  confirms  their1  work;  the 
frequencies  arc  1574*,  175(2),  441(2),  570, 1030, 1357(2),  1394, 
1561(2),  3109  and  3226(2)  cm"1.  For  an  energy  barrier  at  0  K 
(i.e.,  including  zero-point  energy,  ZPE)  of  £0,  the  CTST 
abstraction  rate  coefficient  is 

*(590-1380  K)  =  1.43  x  lO'^T/K)330  x 

exp(-940  K/T)  exp(-E0/RT)  cm3  molecule'1  s'1  (4) 

Preliminary  estimates  showed  that  quantum  mechanical  H-atom 
tunneling  has  little  effect  when  *  >  1  x  1  O'14  cm3  molecule'1  s'1 . 
Agreement  between  eqs  2  and  4  at  825  K,  the  midpoint  of  the 
range  of  T~'  investigated,  is  obtained  when  £0(4)  =  1 2.4  IcJ  mol'1 . 
This  yields  *(  1 380  K)  within  5%  of  the  experimental  value,  while 
*(590  K)  is  too  small  by  a  factor  of  2.3.  Though  accord  at  low 
temperatures  might  be  improved  by  adjustment  of  the  TS 
frequencies,  the  fitted  £0(4)  is  much  lower  than  the  £0(1)  -  67 
kJ  mol'1  from  the  earlier  ab  initio  estimate.1  Use  of  that  value 
together  with  the  observed  total  AlO  removal  rate  coefficient,  eq 
2,  implies  a  branching  ratio  for  direct  abstraction  of  about  9  x 
1 0*3  at  1 380  K  and  about  6X  ICH’at  590  K.  In  order  to  estimate 
the  uncertainty  in  the  £0(1)  value  that  Barve  and  Pettersson 
derived  at  the  CCI  +  Q  level  of  theory,  we  compare  their  estimates 
of  A H°  for  reaction  l  and  the  Li  and  Mg  analogues1  with  the 
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experimental  values.5-6  The  mean  error  is  about  ±25  kJ  mol*1, 
which  is  a  minimum  estimate  of  the  uncertainty  in  £o(l)-  We 
also  allow  for  a  possible  factor  of  3  uncertainty  in  the  TS  partition 
functions.  These  uncertainties  imply  that  branching  ratios  for 
direct  abstraction,  i.e.,  without  formation  of  an  intermediate  bound 
complex,  as  large  as  0.23  at  1380  K  and  0.003  at  590  K  cannot 
be  excluded,  but  another  reaction  channel,  which  is  totally 
dominant  at  the  lower  temperatures,  is  definitely  indicated. 

The  alternative  abstraction  channel  leading  to  HAIO  +  CH3 
would  be  192  ±  93  kJ  mol*1  endothermic5  6  and  can  on  that  basis 
be  excluded  from  further  consideration.  The  likelier  prospects 
for  further  channels  involve  additions  of  AIO  to  CH*.  Particularly, 
insertion  of  AIO  (a  doublet  species)  into  a  C— H  bond  offers  an 
intriguing  prospect,  by  analogy  to  the  suggested  insertion  of  Al 
and  other  group  1 3  doublet  ground-state  atoms  into  such  bonds.22'23 
To  test  this,  we  have  employed  theoretical  methods  to  characterize 
possible  insertion  products  and  estimated  the  likely  pressure 
dependence  of  these  to  see  whether  they  are  consistent  with  our 
observations.  Ab  initio  geometries  of  three  potential  insertion 
species  obtained  by  means  of  MP2/3-21GH  theory  are  shown 
in  Figure  2b-d.  Structure  d  was  previously  proposed  as  the 
product  of  insertion  of  Al  into  CH3OH,24  similar  to  “uch  insertion 
into  H20.24  25  To  confirm  that  these  structures  are  true  energy 
minima,  we  calculated  the  vibrational  frequencies,  which  are  all 
real.26  These  frequencies  are  also  used  to  derive  the  zero-point 
energy  (ZPE),  to  obtain  the  temperature  dependence  of  the 
enthalpies  in  the  kinetic  calculations  outlined  below.  The 
enthalpies  of  CH3AIOH,  CH3A1(0)H,  and  CH3OAlH,  relative 
to  AIO  +  CH4  at  0  K,  are  predicted  at  the  MP4SDQ/6- 
3lG*+ZPE//MP2/3-2lG(*)  level  to  be  -166,  -46,  and  -87  kJ 
mof 1 ,  respectively.27  Equilibrium  constant  calculations  with  these 
enthalpies  indicate  that  CH3A10H  is  stable  under  our  experi¬ 
mental  conditions,  that  the  second  adduct  is  too  weakly  bound 
to  be  a  significant  final  sink  for  AIO,  and  that  formation  of  the 
third  adduct  cannot  be  ruled  out  at  the  lowest  temperatures 
studied.  There  is  a  possibility  that  these  adducts  decompose  to 
new  products.  For  example,  CH3AIOH  might  dissociate  to  AlOH 
+  CH3,  perhaps  with  a  lower  overall  energy  barrier  than  proposed 
for  direct  abstraction. 

Next,  the  likely  pressure  dependence  of  AIO  +  CH4 
CH3AIOH  at  825  K  was  investigated  by  means  of  QRRK  theory.28 
The  addition  rate  coefficient  is  predicted  to  be  within  25%  of  the 
high-pressure  limit  at  1 5  mbar  of  Ar  and  would  show  only  a  small 
pressure  variation  over  the  experimental  range  with  assumed 
Arrhenius  parameters  for  insertion  of  A  —  10~12  cm3  molecule*1 
s- 1  and  E,  =  0.  This  is  consistent  with  the  lack  of  an  observed 
pressure  dependence  of  the  experimental  rate  coefficients.  If 
there  were  a  significant  barrier  to  insertion  (£a  >  0),  then 
dissociation  of  the  excited  adduct  back  to  AIO  +  CH4  would  be 
less  competitive  with  collisional  stabilization  and  the  predicted 
rate  coefficients  would  be  even  closer  to  the  high-pressure  limit. 
The  observed  pressure  independence  is  also  consistent  with  a 
mechanism  where  the  adduct  fragments  to  new  products  before 
stabilization.  Further  QRRK  calculations  suggest  that  formation 
of  any  adduct,  which  is  sufficiently  thermodynamically  stable  to 
be  significant  in  our  experiments,  will  have  a  rate  coefficient  that 
is  in  the  falloff  region  or  close  to  the  high-pressure  limit. 

Conclusions 

These  experiments  arc  in  accord  with  the  suggestion  by  Borve 
and  Pettersson1  that  gas-phase  reactions  between  AIO  and  CH4 
occur  and  may  provide  insight  in  the  heterogeneous  catalytic 
conversion  of  methane  to  C 2  hydrocarbons.  However,  use  of 
their  calculated  barrier  suggests  that  the  branching  ratio  for  direct 
abstraction  is  small  and  therefore  that  other  reaction  channels 
are  important.  Our  calculations  demonstrate  that  insertion 
leading  to  CH3A10H  is  a  reasonable  candidate  for  the  initial 
step,  although  other  channels  are  not  excluded.  The  insertion 


products  (Figure  2b-d)  might  model  the  products  of  dissociative 
adsorption  of  CH4  onto  an  alumina  catalyst,  where  their  energies 
are  likely  to  differ  from  gas-phase  values.  More  detailed 
experimental  and  theoretical  investigations  of  the  present  reaction 
system,  as  well  as  of  other  metal  oxide-hydrocarbon  reactions, 
are  needed  to  obtain  an  understanding  of  this  new  class  of  gas- 
phase  reactions  and  their  potential  importance  for  methane 
conversion. 
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The  high-temperature  fast-flow  reach* ,To  w’”  Tx ' “llTd ifn A 

and  force  constants  that  describe  ,  -  -t  u  obtained  based  on  the  assumption  of  an  0 

N20  reaction  good  agreement  tetwwnMTST  and  ex^n  f^d  for  *e  C02  reaction,  which  indicates 
atom  abstraction  reaction  leading  to  OA1CI.  No  such  agreement  is  i  measured  AlCl  +  02  reaction 

adduct  formation  as  the  mam  AlCl  consumptionchann  .  P  about  150o  K;  however,  adduct 

MTST  calculations  suggest  that  abstraction  can  h*  °* 
formation  appears  to  dominate  over  most  of  the  490-1750  K  range. 


Introduction 

As  part  of  an  ongoing  study  of  the  temperature  dependence 
of  the  rate  coefficients  of  reactions  of  small  Al-  and  B-containmg 
species,  we  report  here  experiments  on  the  reactions 


AlCl  +  N20  —  products 
AlCl  +  C02  —  products 
AlCl  +  SO,  —  products 


0) 

(2) 

(3) 


energy  surfaces,  typical  of  oxidation  reactions  of  metal  specie, 
are  included.  MTST  is  thus  useful  for  such  reactions  when  the 
input  data  for  the  widely  used  semiempirical  1ST  method 
developed  by  Benson,10  principally  for  H/C/O/N  specie  reac¬ 
tions,  are  not  available. 

Experimental  Section 

Detailed  decriptions  of  the  HTFFR  used  in  this  work  have 
been  given  elsewhere.5  6  The  reactor  consists  of  a  vertical  mullite 
reaction  tube  (2.2-cm  i.d.)  surrounded  by  SiC  heating  elements. 
This  assembly  is  sealed  inside  an  insulated  water-cooled  vacuum 
housing.  Two  methods  were  used  to  produce  AlCl.  In  the  first, 
previously  used  method,  a  trace  amount  of  Cl2  or  HC1  was  added 
to  the  Ar  bath  gas  flowing  over  the  Al-wetted  resistively-heated 
tungsten  vaporization  coil,  located  7-17  cm  upstream  fromthe 
reaction  zone.1-3  In  the  second  method,  similar  to  that  used  m 
the  BC1  studies,611  A1C13  was  vaporized  from  a  quartz  crumble 
(4-mm  i  d.)  into  a  stream  of  Ar.  This  stream  then  flowed  through 
a  2450-MHz,  nominally  100-W,  microwave  discharge  in  a  quartz 
tube  ( 1.2-cm  i.d.)  and  then  into  the  reaction  tube.  Both  methods 
were  used  in  the  AlCl  +  C02  study  in  the  1200-1690  K 


We  previously  studied  the  AlCl  reactions  with  02,'  CO:, 2  CV 
and  HC1,4  the  results  of  which  have  been  summarized.  The 
measurements  of  reaction  2  are  repeated  here,  as  the  high- 
temperature  fast-flow  reactor  (HTFFR)  technique  used  has  been 
much  improved5-6  since  the  previous2  study,  where  thedata  showed 
wide  scatter.  The  02  and  C02  reactions  were  originally  thought 
to  lead  directly  to  OA1C1.  However,  a  recent  ab  lmtio  study  of 

the  thermochemistry7  of  OA1C1  suggests  these  two  channels  to  were  used  in  the  AlCl  +  C02  study  in  the  1200-1690  K 

be  too  endothermic8  for  reaction  along  this  path  with  the  observed  temperature  range  and  found  to  yield  identical  results.  The  second 

rate  coefficients.  .  .  method  would  in  general  be  especially  useful  for  HTFFR  studies 

To  investigate  the  mechanism  of  all  these  AlCl  reactions  with  „„aT  rrv,m  temnerature:  such  temperatures  are  not  accessible  by 
oxygen-containing  oxidants,  we  present  here  a  semiempirical 
transition-state  theory  based  method  for  predicting  the  temper¬ 
ature  dependence  of  rate  coefficients.  In  this  model-based 
transition-state  theory  (MTST)  the  rotational  constants  and 
vibrational  frequencies  of  the  transition  state  are  derived  from 
molecular  models  based  on  the  assumption  of  valence  forces. 

This  approach,  the  limitations  of  which  are  discussed,  should  be 
widely  applicable  to  reactions  with  early  activation  barriers  since 
the  required  input  data,  namely,  geometrical  parameters  and 
force  constants,  of  the  transition  state  can  be  estimated  from 
reactant  properties  in  accord  with  Hammond’s  postulate. 

Reactions  occurring  by  a  crossing  of  neutral  and  ionic  potential 


t  Present  address:  Department  of  Chemical  Engineering,  Massachusetts 
Institute  of  Technology,  Cambridge,  MA  02139. 
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near  room  temperature;  such  temperatures  are  not  accessible  by 
the  first  method  because  the  heat  generated  by  the  vaporization 
coil  heats  the  reaction  tube  to  temperatures  greater  than  450  K.. 
Downstream  from  the  source,  a  mixture  of  Ar  and  the  oxidant 
gas  (N2O  CO2,  or  SO2)  was  introduced  into  the  reaction  tube 
through  a'movable  inlet.  The  molar  flow  rate  of  gas  through  the 
inlet  was  kept  constant  during  a  given  rate  coefficient  determi¬ 
nation,  which  consisted  of  five  or  six  measurements  at  varying 
oxidant  concentrations.  A  small  amount  of  02  was  added  to  the 
inlet  gas  to  rapidly5  remove  any  remaining  Al  atoms  entering  the 
reaction  zone.  Such  atoms  otherwise  may  continue  to  produce 
AlCl  and  hence  interfere  with  the  measurements. 

Relative  AlCl  concentrations  were  monitored  by  laser-induced 
fluorescence  using  a  pulsed  Lambda  Physik  EMG  101  excimer/ 
FI  2002  dye  laser  combined  with  a  KDP  doubling  crystal.  The 


0027-3654  ^93/2097-7222504. 00/0  ©  1993  American  Chemical  Society 
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TABLE  I:  Summary  of  Rate  Coefficient  Measurements  for 
A1C1  +  N20*  _  _ 


reaction 


zone 

length 

(cm) 

P 

(mbar) 

[M] 

(1017 

cm-3) 

[N20]m 

(1015 

cm-3) 

F 

(m  s-1) 

T 

(K) 

k  ±  tr*  (cm3 
molecule-1  s-1) 

20  • 

50.9 

4.0 

16.4 

40* 

20 

919 

2.01  ±  0.13  (— 1 4)c 

10 

50.9 

4.0 

16.5 

73 

20 

919 

2.16  ±0.17  (-14) 

20 

50.9 

4.0 

13.7 

36 

21 

921 

2.01  ±0.12  (-14) 

10 

50.9 

4.0 

13.6 

27 

21 

924 

2.18  ±0.19  (-14) 

10 

101.0 

7.7 

14.1 

34 

14 

948 

2.29  ±  0.22  (-14) 

10 

101.2 

7.6 

13.9 

37 

14 

964 

2.08  ±0.22  (-14) 

20 

29.6 

2.2 

14.3 

50 

39 

984 

2.75  ±0.21  (-14) 

20 

29.7 

2.2 

14.2 

62 

39 

991 

2.82  ±0.25  (-14) 

20 

54.1 

5.3 

23.9 

23 

16 

735 

2.40  ±0.33  (-15) 

20 

54.1 

5.4 

24.2 

22 

16 

726 

2.19  ±0.32  (-15) 

20 

78.6 

8.1 

52.3 

30 

7 

704 

1.44  ±0.14  (-15) 

20 

78.6 

8.1 

52.3 

17 

7 

704 

1.80  ±0.13  (-15) 

20 

66.7 

6.9 

40.0 

49 

14 

704 

1.31  ±0.22  (-15) 

20 

66.7 

6.8 

39.7 

48 

14 

709 

1.11  ±0.22  (-15) 

20 

41.6 

4.2 

22.2 

35 

13 

726 

2.23  ±0.22  (-15) 

20 

41.6 

4.1 

22.1 

33 

13 

731 

2.45  ±0.21  (-15) 

20 

55.2 

5.4 

28.9 

34 

10 

740 

2.80  ±0.21  (-15) 

20 

55.2 

5.3 

28.6 

41 

10 

748 

3.13  ±0.22  (-15) 

20 

78.5 

7.2 

35.6 

31 

12 

795 

4.17  ±0.35  (-15) 

20 

62.9 

5.5 

32.7 

50 

17 

830 

7.68  ±0.72  (-15) 

20 

62.9 

5.4 

32.2 

54 

17 

842 

9.51  ±0.64  (-15) 

10 

62.9 

5.5 

32.3 

59 

17 

836 

9.35  ±0.67  (-15) 

10 

49.9 

4.2 

15.0 

23 

19 

856 

1.04  ±0.08  (-14) 

20 

49.9 

4.2 

15.1 

26 

19 

853 

9.80  ±0.70  (-15) 

20 

42.5 

3.6 

13.5 

31 

21 

846 

9.20  ±0.79  (-15) 

20 

42.8 

3.6 

13.5 

24 

21 

852 

1.08  ±0.08  (-14) 

20 

44.3 

3.6 

17.3 

42 

32 

879 

1.30  ±0.16  (-14) 

20 

44.3 

3.6 

17.3 

34 

32 

882 

1.48  ±0.10  (-14) 

20 

35.2 

2.9 

13.9 

17 

40 

873 

1.50  ±  0.15  (-14) 

20 

50.7 

4.2 

20.0 

9 

27 

871 

8.86  ±0.82  (-15) 

20 

68.4 

5.7 

17.3 

52 

16 

875 

1.07  ±0.07  (-14) 

a  The  measurements  are  reported  in  the  sequence  in  which  they  were 
obtained.  b  In  arbitrary  units.  c  Should  be  read  as  (2.0 1  ±  0. 1 3)  x  10-14 
cm3  molecule-1  s-1. 


AlCl^II-X1^)  transition12  was  pumped  on  the  (0,0)  band  at 
261.4  nm,  and  the  fluorescence  intensity  was  observed  through 
a  262-nm  (26-nm  fwhm)  filter  by  an  EMI  9813QA  photomul¬ 
tiplier  tube  coupled  to  a  Data  Precision  Analogic  6000/620  100- 
MHz  transient  digitizer.  Rate  coefficient  measurements  were 
made  in  the  stationary  inlet  mode13  with  observed  reaction  zone 
lengths  of  10  or  20  cm.  The  oxidant  concentrations  were 
maintained  in  great  excess  over  the  A1C1  concentrations.  Rate 
coefficients  and  their  uncertainties  were  obtained  from  the  slopes 
of  linear  plots  of  In  [AlCl]reutive  vs  oxidant  concentration  by  a 
weighted  linear  regression.14 

The  gases  used  were  99.998%  Ar  from  the  liquid,  “Precision 
Aquarator”  99.99%  C02,  Cl2  (1%  in  Ar),  and  “Atomic  Absorp¬ 
tion”  99.0%  N20  from  Linde,  “Ultra  High  Purity”  99.99%  N20 
and  “Anhydrous”  99.98%  S02  from  Matheson,  HC1  (5%  in  He) 
from  Spectra,  and  02  (5.0%  in  Ar)  from  Scott. 

Results 

The  AJG  +  NjO  Reaction.  Reaction  1  was  studied  from  700 
to  990  K.  No  measurements  at  higher  temperatures  are  included 
as  thermal  dissociation  of  N20  is  sufficiently  large  to  interfere 
with  rate  coefficient  measurements.11  Below  700  K  the  reaction 
was  too  slow  to  measure.  Table  I  summarizes  the  rate  coefficients 
and  the  temperatures  at  which  they  were  obtained  along  with  the 
experimental  parameters:  reaction  zone  length,  total  pressure  P, 
total  concentration  [M],  maximum  oxidant  concentration 
[N^O]^,  fluorescence  intensity  F  (a  measure  of  [A1C1]  at  the 
start  of  the  reaction  zone),  and  average  velocity  v.  The 
independence  of  the  individual  rate  coefficients  k{  of  these 
parameters  was  checked  by  examining  plots  of  [fc(  T)  -  k,]/k(  D> 
where  k(T)  is  obtained  from  the  fitted  rate  coefficient  expression 
given  below,  vs  the  above  experimental  parameters.  It  can  be 
seen  from  Figure  1  that  over  the  limited  temperature  range  studied 


r  (K) 

laoo  1000  700 


Figure  1.  Arrhenius  plots  of  experimental  rate  coefficients  for  the  A1CI 
+  N20  (open  circles)  and  A1C1  +  CO2  (open  squares)  reactions.  The 
experimental  rate  coefficient  expressions  are  shown  as  full  lines. 

the  ki  measurements  show  no  deviation  from  Arrhenius  behavior. 
A  weighted  nonlinear  fit15  of  the  data  to  the  expression  k(T)  = 
A  exp(-E/RT)  gives 

*, (700-990  K)  =  5.6  X  10'11  X 

exp(-7380  K/T)  cm3  molecule'1  s'1  (4) 
with  covariance  matrix  elements16 

(<ta/A)2  =  0.1 13,  aE  =  79  900  K2,  o^/A  *  94.6  K 

(5) 

Combining  these  matrix  elements  by  the  propagation  of  errors 
technique410  yields  ±2<r*  precision  limits  of  15%  at  700  K,  6% 
at  850  K,  and  1 2%  at  990  K.  Allowing  10%  uncertainty  associated 
with  the  flow  profile  factor  and  20%  for  further  systematic  errors 
yields  ±2ak  accuracy  limits  of  27%  at  700  K,  23%  at  850,  and 
25%  at  990  K. 

The  AJG  +  COj  Reaction.  Reaction  2  was  studied  over  the 
900-1790  K  temperature  range.  The  rate  coefficient  data 
summarized  in  Table  II  are  again  independent  of  reaction  zone 
length,  P,  [M],  [C02]m«,  F,  and  v.  The  Arrhenius  plot  of  the 
rate  coefficients  shown  in  Figure  1  exhibits  slight  upward 
curvature.  The  rate  coefficient  data  were  fitted  to  the  three- 
parameter  expression  k(T)  =  AT*  exp (-E/RT).  The  goodness 
of  the  fit  x2  is  found  to  be  relatively  insensitive  to  the  value  of 
the  parameter  n  over  the  approximate  range  2.CM.0.  We  therefore 
fix  n  at  3.0.  A  two-parameter  nonlinear  regression  analysis  of 
the  data  gives 

*.,(900-1790  K)  =  4.4  X  10'23  (T/ K)30  X 

exp(-3900  K/T)  cm3  molecule'1  s'1  (6) 
with  covariance  matrix  elements 

(ca/A)2  =  0.0182,  aE2  =  38  000  K2,  a^/A*  26.0  K 

(7) 

Combining  the  covariance  matrix  elements  by  the  propagation 
of  errors  technique  yields  ±2ak  precision  limits  of  17%  at  900  K, 
5%  at  1350  K,  and  6%  at  1790  K.  Allowing  10%  uncertainty 
associated  with  the  flow  profile  factor  and  a  conservative  20%  for 
further  systematic  errors  yields  ±2ak  confidence  limits  of  28% 
at  900  K,  23%  at  1350,  and  23%  at  1790  K.  The  present  rate 
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T4BLE  II:  Summarv  of  Rate  Coefficient  Measurements  for 

Aid  +  ccv  - _ 


reaction 


zone 

length 

(cm) 

P 

(mbar) 

cm-3) 

[co^ 

(1015 

cm-3)  F 

V 

(m 

r>) 

T 

(K) 

k  ±  at  (cm3 
molecule-1  s-1) 

20 

43.3 

2.5 

18 

18* 

25 

1273 

4.20 

±0.55  (-15)^ 

20 

78.4 

4.5 

33 

71 

13 

1254 

3.73 

±0.50  (-15)' 

20 

78.6 

4.5 

32 

64 

14 

1273 

3.21 

±0.22  (-15)' 

10 

78.6 

4.4 

32 

62 

14 

1281 

3.86 

±0.53  (-15)' 

10 

78.6 

4.4 

32 

108 

14 

1284 

4.15 

±0.2*  (-15)' 

10 

63.5 

3.2 

17 

33 

14 

1429 

8.71 

±0.9*  (-15)' 

10 

63.7 

3.2 

21 

47 

14 

1440 

9.99 

±0.91  (-15)' 

20 

62.8 

3.1 

12 

50 

14 

1475 

9.69 

±0.9*  (-15)' 

20 

62.8 

3.0 

12 

54 

14 

1493 

1.28 

±  0.13  (-14)' 

20 

43.5 

2.8 

16 

13 

19 

1132 

2.84 

±0.30  (-15)' 

20 

43.5 

2.7 

15 

18 

19 

1149 

2.38 

±0.9*  (-15)' 

20 

58.1 

4.3 

23 

10 

13 

990 

7.83 

±  2.20  (-16)' 

20 

57.7 

4.1 

23 

23 

13 

1010 

1.06 

±0.1*  (-15)' 

20 

88.1 

6.2 

27 

27 

11 

1031 

1.15 

±0.17  (-15)' 

20 

46.4 

3.1 

20 

59 

22 

1091 

1.27 

±0.34  (-15)' 

20 

46.1 

3.1 

20 

48 

22 

1094 

1.34 

±0.45  (-15)' 

20 

53.5 

2.4 

14 

21 

31 

1584 

1.33 

±0.19  (-14)' 

10 

53.3 

2.4 

14 

24 

31 

1602 

1.74 

±0.21  (-14)' 

10 

53.2 

2.4 

14 

21 

31 

1609 

2.04 

±0.18  (-14)' 

10 

92.0 

4.0 

20 

21 

15 

1679 

i.tv 

±0.24  (-14)' 

10 

92.0 

4.0 

20 

17 

15 

1686 

2.05 

±0.21  (-14)' 

20 

53.7 

4.2 

21 

55 

11 

930 

7.20 

±  3.3*  (-16)' 

20 

53.7 

4.1 

21 

43 

11 

942 

9.32 

±3.11  (-16)' 

20 

65.3 

5.1 

27 

74 

16 

920 

5.70 

±  2.62  (-16)' 

20 

65.7 

5.1 

27 

54 

16 

925 

6.69 

±  1.47  (-16)' 

20 

78.6 

6.3 

36 

39 

16 

900 

5.97 

±  1.34  (-16)' 

20 

43.1 

2.3 

17 

19 

34 

1368 

5.49 

±  1.3*  (-15)' 

20 

43.1 

2.3 

17 

40 

35 

1373 

5.24 

±0.71  (-15)' 

10 

43.1 

2.3 

17 

68 

35 

1374 

6.10 

±  2.76  (-15)' 

10 

43.1 

2.3 

17 

66 

35 

1386 

6.77 

±0.83  (-15)' 

10 

85.8 

4.5 

33 

40 

18 

1397 

6.66 

±0.6*  (-15)' 

10 

85.4 

4.4 

33 

19 

18 

1399 

6.00 

±0.95  (-15)' 

20 

83.6 

4.3  - 

17 

30 

18 

1401 

5.74 

±  1.00  (-15)' 

20 

99.6 

6.0 

30 

113 

18 

1210 

3.69 

±0.41  (-15)' 

20 

113.3 

6.8 

30 

87 

18 

1204 

2.72 

±0.2*  (-15)' 

20 

73.3 

4.2 

25 

59 

17 

1250 

3.75 

±0.32  (-15)' 

20 

60.4 

2.8 

12 

36 

45 

1540 

1.81 

±0.15  (-14)' 

20 

60.4 

2.8 

12 

39 

46 

1556 

1.80 

±0.23  (-14)' 

20 

91.4 

4.2 

16 

39 

25 

1590 

1.83 

±0.10  (-14)' 

20 

91.6 

4.2 

16 

78 

25 

1590 

1.90 

±0.11  (-14)' 

10 

92.8 

4.2 

21 

53 

25 

1599 

1.71 

±0.12  (-14)' 

10 

93.0 

4.2 

21 

54 

25 

1599 

1.74 

±0.12  (-14)' 

10 

61.9 

2.7 

18 

62 

22 

1632 

2.14 

±0.15  (-14)' 

20 

60.5 

2.7 

12 

65 

22 

1638 

2.32 

±  0.14  (-14)' 

20 

52.9 

2.6 

17 

87 

33 

1448 

6.59 

±0.54  (-15)' 

20 

52.8 

2.6 

16 

59 

33 

1449 

6.98 

±0.64  (-15)' 

10 

52.8 

2.6 

16 

60 

33 

1454 

9.33 

±1.11  (-15)' 

10 

52.8 

2.6 

16 

60 

33 

1454 

9.08 

±  1.06  (-15)* 

10 

75.0 

3.1 

12 

70 

37 

1765 

2.75 

±0.28  (-14)' 

10 

75.2 

3.1 

12 

69 

36 

1765 

2.81 

±  0.20  (-14)' 

20 

74.8 

3.0 

7 

87 

37 

1777 

2.26 

±0.14  (-14)' 

20 

74.8 

3.0 

12 

88 

37 

1777 

2.57 

±0.24  (-14)' 

20 

41.9 

1.7 

8 

91 

68 

1787 

2.61 

±0.21  (-14)' 

20 

42.1 

1.7 

8 

125 

67 

1787 

2.55 

±0.16  (-14)' 

20 

42.1 

1.7 

8 

39 

67 

1791 

2.98 

±0.26  (-14)' 

«  The  measurements  are  reported  in  the  sequence  in  which  they  were 
obtained.  b  In  arbitrary  units.  *  Should  be  read  as  (4.20  ±  0.55)  x  10-15 
cm3  molecule-1  s-1 .  4  A1C1  produced  by  the  microwave-discharge  method. 
f  .4101  produced  by  the  reacting  CI2  or  HC1  with  an  Al-wetted  vaporization 

coil. 


coefficient  recommendation  is  in  good  agreement  with  that 
determined  in  the  previous  investigation;2  however,  the  scatter  in 
the  measurements  is  now  greatly  reduced. 

The  A1G  +  SO*  Reaction.  The  rate  coefficients  of  reaction 
3  could  not  be  accurately  determined  because  of  absorption  and 
emission  of  laser  radiation  by  S02,12  which  greatly  lowered  the 
sensitivity  of  A1C1  detection.  Corrections  were  made  to  the 
integrated  fluorescence  signal  to  remove  the  S02  contribution. 
Over  the  800-1 100  K  temperature  range  no  significant  change 
in  the  A1C1  signal  could  be  observed  as  the  S02  concentration 
was  varied,  which  indicates  an  upper  limit  to  the  rate  coefficient 


Jt3(800— 1 100  K)  <  5  X  10_M  cm3  molecule-1  s-1  (8) 

We  were  unable  to  investigate  reaction  3  at  higher  temperatures 
because  the  emitted  radiation  from  species  other  than  A1C1  was 
too  intense.  Thermochemical8  equilibrium  calculations  indicate 
that  above  1 100  K  SO,  which  also  absorbs  and  emits  radiation 
near  the  A1C1  transition,12  can  be  present  in  significant  quantities. 

Model-Based  Transition-State  Theory  (MTST) 

According  to  conventional  transition-state  theory  the  rate 
coefficient  for  a  reaction  occurring  on  a  single  potential  energy 
surface  is  given  by17 

k(T)  *  nr  TT  «P  RT)  (9) 

in  which  q\  is  the  partition  function  for  the  transition  state,  qA 
and  qs  are  partition  functions  for  the  reactants,  and  A£*  is  the 
difference  in  zero-point  levels  of  the  transition  state  and  the 
reactants.  qA  and  ?b  are  readily  evaluated  from  experimental 
values  of  the  rotational  constants  and  vibrational  frequencies  of 
the  reactants.  The  data  needed  to  evaluate  q\  are  seldom  available 
and  must  therefore  be  estimated.  The  MTST  approach  for  this 
evaluation  is  suitable  for  transition  states  which  arise  from  an 
early  crossing  of  neutral  and  ionic  potential  surfaces,  such  as  is 
the  case  for  many  metal  species  reactions.  For  such  a  transition 
state  reasonably  good  estimates  of  the  geometrical  parameters 
and  vibrational  frequencies  can  be  made  using  Hammond’s 
postulate.9  It  can  then  be  assumed  that  the  Intemuclear  distances, 
valence  angles,  and  force  constants  of  the  transition  states  are 
closely  related  to  the  properties  of  the  reactants  for  an  exothermic 
reaction  and  properties  of  the  products  for  an  endothermic 
reaction. 

Vibrational  frequencies  for  transition  states  are  obtained  from 
molecular  models.  The  molecular  model  used  here  assumes  that 
harmonic  restoring  forces  oppose  changes  in  bond  lengths,  angles 
between  valence  bonds,  and  dihedral  angles.  The  potential  energy 
V  in  the  region  near  the  transition  state  is  thus  given  by18 

2V=  JVrfAr,*  +  A*/  +  X/,*Ar*2  (10) 

i  J  * 

where  Fn  F*,  and  Fr  are  force  constants  for  changes  in  intemuclear 
distances  r,  valence  angles  6,  and  dihedral  angles  a  from  their 
respective  equilibrium  values.  A  normal-coordinate  analysis  of 
the  model,  which  procedure  has  been  described  in  detail  by  Pauling 
and  Wilson,19  is  carried  out.  The  secular  equation  thus  obtained 
is  solved  numerically  by  reducing  the  determinant  to  the 
tridiagonal  form  using  Householder’s  algorithm.15  The  char¬ 
acteristic  polynomial  of  the  resulting  matrix  is  then  evaluated  by 
the  QL  algorithm  with  implicit  shifts.15  Vibrational  frequencies 
Vi  are  obtained  from  the  relation19 

= x^/itc  do 

in  which  A,  are  the  eigenvalues  of  the  secular  equation .  Rotational 
constants  Bt  for  the  transition  state  are  obtained  from  the 
eigenvalues  of  the  inertia  tensor,20  that  is,  from  the  principal 
moments  of  inertia  /,  by  the  relation21 

5,  =  h/C&^cI,)  (12) 

Several  examples  of  accordingly  estimating  transition-state 
parameters  for  A1C1  reactions  are  given  in  the  Discussion. 

Some  limitations  of  MTST  are  that  force  constants  arising 
from  the  reaction,  which  includes  one  stretching  force  constant 
and  one  or  more  bending  or  torsional  force  constants,  must  be 
estimated  from  the  available  ab  initio  calculations  on  similar 
reactions.  The  most  important  of  these  are  the  torsional  force 
constants,  since  the  torsions  can  be  of  low  frequency  and  hence 
influence  the  TST  rate  expression.  Also,  MTST  is  not  suitable 
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TABLE  III:  Transition-State  Parameters  for  C1A10NN* 

geometrical  parameter^  force  constant* 


r(ClAi) 

2.13 

FV(CIAI) 

2.1 

r(AlO) 

2.25 

Fri  AlO) 

-2.0 

r(ON) 

1.18 

Fr(  ON) 

5.85 

r{  NN) 

1.13 

Fri  NN) 

18.7 

0(C1A1O) 

120 

Fj(ClAJO) 

0.2 

0(A1ON) 

90 

Ff(AlON) 

0.6 

0(ONN) 

160 

Fj(ONN) 

0.67 

r(CLAlON) 

0 

FT(C1A10N) 

0.02 

r(  AlONN) 

0 

Fri  AlONN) 

0.02 

rotational  constants  (B \BjBi):  9.51  x  10-4  cm-3 
vibrational  frequencies  y,:  39,  89,  150,  291,  478,  605,  991,  2221  cm-1 

We  assume  trans  arrangements  for  the  atoms  in  the  C1A10N  and 
AlONN  groups. b  Units  of  r  are  10~10  m;  units  of  6  and  r  are  deg. c  Units 
of  Fr  are  102  N  nr1;  units  of  F 9  and  Fr  are  10“18  N  m  rad*2. 

for  reactions  with  a  centrally  located  transition  state,  e.g.,  H  + 
HC1  —  H2  +  Cl,  because  the  potential  energy  function  (cf.  eq 
10)  would  have  to  be  modified  to  include  a  cross- term  for  the 
interaction  between  the  H-H  and  H-Cl  bonds.22  Furthermore, 
the  transition  state  would  resemble  both  the  rear  ints  and  the 
products  in  part;  hence,  there  would  be  considerable  uncertainty, 
in  estimating  geometrical  parameters  and  force  constants. 

Discussion 

The  A1C1  +  N20  Reaction.  The  lowest  channel  of  the 
abstraction  reaction 

A1C1  +  N.O  —  O  A1C1  +  N2  (la) 

is  spin-allowed  and  267  kJ  mol-1  exothermic.23  The  parameters 
of  the  transition  state  CIAIONN*  are  derived  from  properties  of 
the  reactants  and  are  summarized  in  Table  III.  The  distances 
r(ClAl),  r( ON),  and  r(NN)  are  from  values  for  the  A1C1  and 
N20  molecules.  The  intemuclear  distance  along  the  reaction 
coordinate  r(A10)  is  estimated  by 

r(AlO)  =  [r(Al+)  +  r(CT)]0.50,/(i-1)  (13) 

Equation  13  is  based  on  a  relation  given  by  Pauling24  between 
univalent  radii  and  the  bond  distances  of  diatomic  alkali-metal 
halide  molecules.  We  assume  that  the  relation  applies  to  charge- 
transfer  complexes  and  hence  to  transition  states,  which  are  located 
fairly  close  to  the  complexes.  Values  for  the  univalent  radii  r{  Al+) 
and  /•(<>)  and  the  average  Born  exponent  5  are  obtained  as 
described  in  ref  25.  The  constant  0.50  is  obtained  such  that  eq 
1 3  gives  good  agreement  with  ab  initio  computations26-2*  of  B-O, 
Al-O,  and  H-O  distances  for  transition  states  of  the  respective 
BC1  +  02,  A1  +  C02,  and  H  +  N20  reactions.  Lacking 
information  on  the  valence  and  dihedral  angles  of  the  transition 
state,  we  assume  the  values  given  in  Table  III  based  on  optimized 
geometries  from  ab  initio  computations.26-28  The  force  constant 
Fr(C IA1)  is  taken  as  that  in  A1C1,29  while  Fr(NN)  and  F*(ONN) 
are  obtained  from  data  for  N2O.30  The  remaining  force  constants 
in  Table  III  are  estimated  by  a  scheme  that  gives  good  agreement 
with  ab  initio  computations26-28  of  vibrational  frequencies.  For 
example,  Fr(ON)  is  taken  as  one-half  the  N-O  stretching  force 
constant  in  N20.  In  addition,  all  torsional  force  constants  are 
taken  as  0.02  x  1 0~18  N  m  rad-2  and  the  stretching  force  constants 
corresponding  to  reaction  coordinates  are  taken  as  -2.0  x  102  N 
nr1.  These  data  are  used  to  calculate  the  transition-state 
rotational  constants  and  vibrational  frequencies,  which  are 
summarized  in  Table  III. 

To  allow  comparison  between  theory  and  experiment,  the 
preexponential  factor  of  eq  9  is  calculated  for  the  500-2000  K 
temperature  range  and  fitted  to  the  expression  A  T \  To  calculate 
the  barrier  height,  the  experimental  rate  coefficient  expression 
is  set  equal  to  the  TST  equation  (9)  at  820  K,  the  temperature 
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Figure  2.  Arrhenius  plots  of  experimental  rate  coefficients  for  the  A1CI 
+  N20  (open  circles)  and  A1C1  +  C02  (open  squares)  reactions.  Rate 
coefficient  expressions  for  the  O  atom  abstraction  channels  for  these 
reactions  derived  using  the  MTST  method  are  shown  as  solid  lines. 

TABLE  IV:  Transition-State  Parameters  for  CIA  IOC  O* 

geometrical  parameter^  _ force  constant* _ 


KClAl) 

2.08 

FAC1A1) 

3.2 

KAIO) 

1.63 

F^AIO) 

4.2 

HOC) 

1.81 

Fr(OC) 

-2.0 

r(C  O) 

1.13 

FriC  O) 

19.0 

0(C1A1O) 

160 

F,(CIA10) 

0.3 

0(A1OC) 

90 

F#(AiOC) 

0.6 

6(OCO) 

120 

F,(OCO) 

0.2 

r(ClAlOC) 

0 

FT(CIA10C) 

0.02 

r(  AlOCO) 

0 

Fr(  AlOCO) 

0.02 

rotational  constant s  C 5i52B3):  9.46  x  10-4  cm-3 
vibrational  frequencies  y,-:  75,  91,  141,  216,  231,  523,  908,  2164  cm-1 

a  WTe  assume  trans  arrangements  for  the  atoms  in  the  CIAIOC  and 
AlOCO  groups.  The  parameters  of  the  transition  state  are  based  on  the 
products  of  the  endothermic  reaction  A1C1  +  C02  —  CO  +  OA1C1. 
b  Units  of  r  are  10-10  m:  units  of  6  and  r  are  deg.  c  Units  of  Fr  are  102 
N  m-1;  units  of  F$  and  Fr  arc  10-18  N  m  rad-2. 

corresponding  to  the  middle  of  the  T~l  range  of  the  experiments. 
The  value  A E1  =  45.9  kJ  mol-1  is  thus  obtained,  which  yields 

*MTST  =  5  74  X  10‘18  {T/K)1M  X 

exp(-5520  K/7")  cm3  molecule'1  s'1  (14) 

A  comparison  between  eq  14,  the  MTST  rate  coefficient  for 
reaction  la,  and  experiment  is  shown  in  Figure  2.  Over  the 
temperature  range  of  the  experimental  data  the  agreement  can 
be  seen  to  be  excellent,  which  suggests  that  the  dominant  channel 
for  reaction  1  is  abstraction,  i.e.,  channel  la. 

The  A1C1  +  CO2  Reaction.  The  abstraction  reaction 

A1C1  +  C02  —  OA1C1  +  CO  (2a) 

is  98  kJ  mol-1  endothermic.23  The  parameters  of  the  transition 
state  CIAIOCO2,  which  are  summarized  in  Table  IV,  have  been 
derived  based  on  properties  of  the  products  OA1C1  +  CO.  The 
estimation  scheme  is  very  similar  to  that  for  the  A1C1  +  N20 
transition  state  so  only  the  differences  will  be  discussed  here.  The 
r(CO),  KC1A1),  and  r( AlO)  distances  are  taken  as  those  in  CO 
and  OA1C1.  The  intemuclear  distance  along  the  reaction 
coordinate  is  estimated  by  r(OC)  =  [KO“)  +  r(C+)]0.501/<5-1>. 
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The  force  constant  F,(CO)  is  taken  as  that  in  the  CO  molecule, 
while  F,(CLA1)  and  F»(C1A10)  are  taken  as  equal  to  force  constants 
derived  from  ab  initio  calculations7  of  the  vibrational  frequencies 
of  O  A1C1.  Table  IV  summarizes  the  resulting  rotational  constants 
and  vibrational  frequencies. 

We  calculate  the  preexponential  factor  of  eq  9  in  the  500- 
2000  K  temperature  range  and  fit  the  results  to  the  expression 
A  V1  expression.  A  barrier  height  of  1 1 8  kJ  mol-1  for  the  forward 
reaction  is  estimated  based  on  the  endothermicity  of  reaction  2a 
plus  a  barrier  height  for  the  reverse  reaction  of  20  kJ  mol'1, 
which  is  probably  too  small  in  view  of  the  fact  that  CO  +  OA1C1 
is  isoelectronic  with  A1C1  +  N20,  for  which  we  obtained  a  barrier 
height  of  45.9  kJ  mol'1 .  Combining  the  calculated  A ,  n,  and  A£* 
values  leads  to  the  expression 

*MTST  =  I'27  X  10"18  (W2'44  X 

exp(_14190  K/7)  cm3  molecule"1  s'1  (15) 

Equation  15  is  plotted  in  Figure  2  to  allow  for  comparison  with 
experiment.  It  may  be  seen  that  the  rate  coefficients  thus 
calculated  are  too  low  to  account  for  the  experimental  data  below 
about  1 400  K.  If  a  higher  barrier  for  the  reverse  ruction  were 
to  be  included,  the  calculated  rate  coefficients  would  decrease 
further,  reinforcing  this  conclusion.  We  thus  speculate  that  at 
least  below  1400  K  the  main  products  of  reaction  2  are  adducts 
correlating  with  the  singlet  and  triplet  potential  energy  surfaces 
of  A1C1+  with  C02" 

The  A1G  +  SO*  Reaction.  The  abstraction  reaction 

A1C1  +  S02  —  OA1C1  +  SO  (3a) 

has  an  endothermicity23  of  119  kJ  mol"1  and  is  spin- forbidden. 

For  temperatures  between  800  and  1 100  K  the  rate  coefficient 
for  abstraction  is  much  lower  than  the  upper  limit  determination 
(cf.  eq  8)  because  of  the  large  endothermicity.  There  may  be 
channels  leading  to  the  formation  of  adducts.  However,  since  no 
reaction  could  be  detected,  the  adducts  must  have  low  binding 
energies  and  the  addition  channels  may  have  significant  activation 
barriers. 

The  A1C1  +  O2  Reaction.  In  the  study1  of  the  rate  coefficients 
of  the  A1C1  +  02  reaction  from  490  to  1 750  K,  it  was  speculated 
that  the  dominant  low  temperature  channel  could  be 

A1C1  +  02  — *  A102  +  Cl  (16a) 

while  at  higher  temperatures  the  channel 

A1C1  +  02  —  OA1C1  +  0  (16b) 

was  considered  a  possibility.  Thermochemical  data  available  at 
that  time  indicated  that  the  thermochemistry  of  channel  1 6a  was 
uncertain,  while  channel  16b  was  thought  to  be  exothermic.1 
Recent  high-level  ab  initio  computations7  on  A102  now  indicate 
that  channel  16a  is  87  kJ  mol"1  endothermic,  confirming  the 
A//KAIO2)  assumed  in  the  JANAF  tables.8  That  theoretical 
work23  also  suggests  that  channel  1 6b  is  66  kJ  mol"1  endothermic. 
These  data  and  the  use  of  MTST  allow  a  reevaluation  of  the 
measurements  of  the  A1C1  +  O2  reaction. 

The  parameters  for  the  molecular  model  as  well  as  the  rotational 
constants  and  vibrational  frequencies  of  the  transition  state  for 
channel  16a  are  summarized  in  Table  V.  The  geometry  and 
force  constants  have  been  obtained  based  on  the  assumption  that 
the  transition  state  resembles  the  products  OAlO  +  Cl,  with  the 
halogen  atom  attacking  the  A1  atom  by  a  side-on  approach.  We 
assume  the  AlO  distance  and  force  constants  Fr(A10)  and  Fr 
(OAIO)  for  the  transition  state  are  equal  to  those  of  A102.8  The 
remaining  parameters  have  been  obtained  similarly  to  those  for 
the  A1C1  +  N20,  C02  reactions.  If  the  barrier  height  for  the 


TABLE  V:  Transition-State  Parameters  for  AlO^CF 

geometrical  parameter0  force  constant6 


r(AlCl) 

2.34 

FJ  A1C1) 

-2.0 

rfAlO) 

1.69 

FJAIO) 

4.2 

0(OA1O) 

160 

FtiOAXO) 

0.24 

fl(OAlCl) 

100 

F^OAICI) 

0.2 

t(A102C1) 

0 

FT(A102CI) 

0.02 

rotational  constants  (S1B2B3):  1.72  x  10"3  cm-3 
vibrational  frequencies  1 19, 176,  200,  678,  981  cm"1 


a  We  assume  a  arrangement  for  the  atoms.  The  parameters  of  the 
transition  state  are  based  on  the  products  of  the  endothermic  reaction 
A1C1  +  O:  —  OAlO  +  Cl.  b  Units  of  r  are  10"10  m;  units  of  0  and  r  are 
deg. c  Units  of  Fr  are  102  N  nr1;  units  of  F#  and  Fr  are  1(HS  N  m  rad-2. 

TABLE  VI;  Transition-State  Parameters  for  CIAIOO* 
geometrical  parameter0  force  constant* 


KC1AI) 

2.08 

F^ClAl) 

3.2 

r(A10) 

1.63 

FlfAlO) 

4.2 

rioo) 

1.72 

FAOO) 

-2.0 

fl(ClAlO) 

160 

/V<C1A10) 

0.3 

fl(AlOO) 

90 

/V(A100) 

0.6 

t(CIAIOO) 

0 

/■.(CIAIOO) 

0.02 

rotational  constants  3.27  X  10"3  cm-3 

vibrational  frequencies  n,:  125, 131, 245, 522,  907  cm"1 

a  We  assume  a  trans  arrangement  for  the  atoms.  The  parameters  of 
the  transition  state  are  based  on  the  products  of  the  endothermic  reaction 
A1C1  +  o2  —  O  +  OA1CL  b  Units  of  r  are  10-10  m;  units  of  0  and  r  are 
deg. e  Units  of  Fr  arc  102  N  nr1;  units  of  F»  and  Fr  are  10~1$  N  m  rad"2. 

T  (K) 


Figure  3.  Arrhenius  plot  of  experimental  rate  coefficients  for  the 
previously1  studied  A1C1  +  02  (open  circles)  reaction.  Rate  coefficient 
expressions  for  the  A1C1  +  02  —  OA1C1  +  O  (full  line)  and  A1C1  +  O* 
A102  +  Cl  (dashed  line)  channels  derived  using  the  MTST  method 
are  shown  for  comparison  with  experiments. 

reverse  reaction  is  neglected,  MTST  calculations  give 

*MTST  “  2-75  X  10-17  (TV K)1-76  X 

exp(- 10460  K/7")  cm3  molecule"' s"1  (17) 

As  can  be  seen  in  Figure  3  the  MTST  calculations  given  by  the 
dashed  line  indicate  that  channel  16a  is  insignificant  at  all 
temperatures. 

The  MTST  parameters  for  channel  16b  are  summarized  in 
Table  VI  and  have  been  obtained  similarly  to  those  for  the  A1C1 


Reactions  of  A1C1  with  N20.  C02,  and  S02 

+  C02  reaction.  If  we  assume  a  negligible  barrier  height  for  the 
reverse  reaction,  the  MTST  expression  obtained  is 

W=U8X10-17(W‘87X 

exp(-7938  K/T)  cm3  molecule'1  s'1  (18) 

In  Figure  3  the  experimental  rate  coefficient  measurements  are 
shown  along  with  those  from  eq  18.  This  comparison  suggests 
that  channel  16b  can  be  of  significance  above  about  1500  K  in 
agreement  with  the  earlier  speculation.1 

Recent  ab  initio  calculations26  have  indicated  the  existence  of 
a  stable  adduct  A1C102,  which  correlates  with  the  reactants  A1C1 
+  o2.  Hence,  the  dominant  reaction  pathway  at  temperatures 
below  1 750  K  is  likely  to  be  addition  with  the  abstraction  reaction 
1 6b  becoming  of  greater  importance  with  increasing  temperature. 
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Rate  coefficients  for  the  title  reactions  have  been  measured  in  a  high-temperature  fast-flow  reactor  (HTFFR). 
Weighted  fits  give  for  BC1  +  S02  *(470-1690  K)  =  7.7  X  lO-^T/K)30  exp(-1680  K/D  cm3  molecule'1  r1 
and  for  BC1  +  N20  *(670-1000  K)  =  5.1  x  10~12  exp(-6340  K/7)  cm3  molecule'1  r>.  2tr  accuracy  limits 
are  about  ±25%.  Model-based  transition-state  theory  (MTST)  is  used  to  calculate  rate  coefficients  for  the 
O-atom  abstraction  channels  of  these  reactions,  as  well  as  for  the  previously  measured  BC1  reactions  with  02 
and  CO2.  Approximate  agreement  between  theory  and  experiment  is  obtained  for  the  four  reactions,  indicating 
abstraction  to  be  dominant  for  each  reaction.  The  reactions  are  compared  to  the  corresponding  A1C1  reactions. 


Introduction 

In  previous  work,  the  high-temperature  fast-flow  reactor 
(HTFFR)  technique  has  been  used  to  make  rate  coefficient 
measurements  on  the  reactions  BC1  +  02,  in  the  540-1670  K 
range,1  and  BC1  +  C02,  over  the  770-1830  K  range.2  Here  we 
extend  this  series  by  reporting  measurements  on  the  reactions 


BC1  +  S02  — ■  products 

(i) 

BC1  +  N20  —  products 

(2) 

in  the  respective  470-1690  K  and  670-1000  K  temperature  ranges. 

In  another  recent  paper,  we  reported  on  the  A1C1  reactions 
with  these  same  four  oxygen-containing  oxidants  and  developed 
a  variant  of  classical  transition-state  theory  (CTST)  .3  This  model- 
based  transition-state  theory  (MTST)  is  applicable  to  reactions 
with  early  activation  barriers  and  is  particularly  useful  when  the 
input  data  for  the  CTST  methods  originated  by  Benson4  are  not 
available.  MTST  was  used  in  that  work3  to  determine  the  AT1 
part  of  the  rate  coefficient  k(T)  =  AT1  exp(-£ / RT)  expressions 
for  abstraction  reactions.  Here  we  continue  the  comparison 
between  experiment  and  MTST  calculations  for  reactions  of  BC1, 
isoelectronic  with  AiCl. 

Experimental  Section 

The  HTFFR  technique,  as  used  for  BC1  studies,  has  been 
described  in  detail  elsewhere,2'5  Briefly,  BC1  was  generated  by 
passing  20-40  ppm  BCl3/Ar  through  a  microwave  discharge  in 
a  Pyrex  tube  (1.2-cm  i.d.)  and  fed  into  the  reaction  tube.  The 
reactor  consisted  of  a  vertically  positioned  mullite  tube  (2.2-cm 
i.d.)  sealed  inside  a  water-cooled,  steel  vacuum  housing.  Two 
parallel  series  of  SiC  heating  elements  were  used  to  maintain 
uniform  temperatures  in  the  reaction  tube.  Efficient  heating  of 
the  bath  gas  was  achieved  by  fixing  the  smaller  mullite  tube 
(1.8-cm  o.d.)  in  the  reaction  tube  upstream  from  the  reaction 
zone.  The  oxidant  gas,  S02  or  N20  mixed  with  N 2,  was  introduced 
into  the  reaction  tube  through  a  movable  alumina  ring  inlet  with 
eight  circular  openings  (0.1-cm  diameter)  through  which  the 
oxidant  gas  mixtures  were  injected  into  the  BC1/  Ar  stream.  Gas 
flow  through  the  inlet  was  maintained  at  a  constant  molar  rate 
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for  rate  coefficient  determinations,  each  of  which  consisted  of 
five  or  six  measurements  at  different  oxidant  concentrations.  At 
10  or  20  cm  downstream  from  the  ring,  relative  concentrations 
of  BC1  were  monitored  by  laser-induced  fluorescence  (LIF)  using 
four  different  transitions  of  the  (A^-X1^)  band,  namely,  ( 1 ,0)Q 
at  266.0  nm,  (0,0)Q  at  272.0  nm,  (0,0)P  at  272. 16  nm,  and  (1 ,1  )P 
at  272.2 1  nm.6  Laser  pulses  at  these  wavelengths  were  generated 
by  a  Lambda  Physik  EMG 101  excimer  laser/ FL  2002  dye  laser/ 
KDP  doubling  crystal  combination.  The  pulses  were  directed 
through  two  4-mm-diameter  holes  in  the  mullite  reaction  tube. 
The  BC1  fluorescence,  which  passed  through  a  third  hole  in  the 
reaction  tube  positioned  normal  to  the  laser  pulse  path,  was 
recorded  with  an  EMI  9813  QA  photomultiplier  tube  (PMT )/ 
LeCroy  VV100B  widc-band  pulse  amplifier /Data  Precision 
Analogic  6000/620  100-MHz  transient  digitizer  combination. 
A  270-nm  interference  filter  (24-nm  fwhm)  was  positioned  in 
front  of  the  PMT  to  remove  background  radiation. 

Rate  coefficient  measurements  were  made  in  the  stationary 
inlet  mode7  with  reaction  zone  lengths  of  10  or  20  cm,  under 
pseudo-first-order  conditions  with  BC1  as  the  minor  reactant. 
Plots  of  In  [BC1]  relative  VS  [S02]  OT  [N20],  in  which  [BC1] relative 
is  the  time- integrated  LIF  signal  averaged  over  200  laser  pulses, 
were  linear  with  slopes  proportional  to  the  rate  coefficients.  For 
each  individual  measurement,  kj  and  a*,  were  determined  by  a 
weighted  linear  least-squares  fit  of  the  data.7’8 

The  gases  used  were  Ar  (99.998%,  Linde)  and  N2  (99.998%, 
Airco)  obtained  from  the  liquid,  BC13  (99.9%,  Matheson),  N20 
(99.99%,  Matheson),  N20  (99.0%,  Linde),  S02  (99.98%,  Mathe¬ 
son),  S02  (99.98%,  Linde),  and  02  (5.0%  in  Ar,  certified;  Scott). 

Results 

The  measured  rate  coefficients  k{  of  reactions  1  and  2  at  the 
various  temperatures  used  are  shown  in  Figure  1.  They  are 
summarized  in  Tables  I  and  II  along  with  the  observed  reaction 
zone  length,  total  pressure  P,  total  concentration  [M],  maximum 
oxidant  concentration  [S02]nn*  or  [N^jn^x,  fluorescence  intensity 
F  (a  measure  of  [BC1]  at  the  start  of  the  reaction  zone),  and 
average  velocity  v.  Variation  in  the  experimental  conditions  was 
used  to  check  the  accuracy  of  the  k,  measurements,  as  secondary 
reactions  of  potential  other  species  produced  in  the  Ar  +  BC1 
discharge,  i.e.,  B,  BC1,  BC12,  Cl,  and  Cl2,  could  have  produced 
BC1  in  the  reaction  zone.  If  this  had  occurred,  k{  measurements 
would  have  shown  dependences  on  reaction  zone  length,  t),  and 
P.  Our  reported  fc,  values  (Tables  I  and  II)  may  be  seen  to  be 
independent  of  these  quantities  by  examining  plots  of  [k(T)  - 
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Figure  1.  Arrhenius  plots  of  experimental  rate  coefficients  for  the  BC1 
+  S02  (open  circles)  and  BC1  +  N20  (open  squares)  reactions.  The  full 
lines  represent  the  experimental  rate  coefficient  expressions. 

ki]/k(T)y  where  k(T)  are  the  fitted  rate  coefficient  expressions 
given  below,  vs  the  experimental  parameters;  for  each  case,  the 
ki  data  are  randomly  scattered  about  the  zero  line  of  the  ordinate. 

BC1  +  S02  Reaction.  The  85  rate  coefficient  measurements 
of  this  reaction  were  obtained  over  the  temperature  range  470- 
1690  K  (Table  I).  In  this  study,  LIF  radiation  was  observed  in 
the  absence  of  BC1  at  temperatures  above  about  1300  K.  For 
each  measurement  with  BCI,  the  intensity  of  this  radiation  was 
subtracted  from  the  signal.  The  interference  was  observed  to 
increase  as  the  reaction  tube  temperature  was  raised  and  as  the 
[S02]  increased.  Since  S02  as  well  as  SO,  generated  from  thermal 
dissociation  of  S02,  absorb  in  the  260-280-nm  region,9  both  could 
have  contributed  to  this  background  radiation.  There  was  also 
an  a  priori  possibility  that  SO,  produced  by  reaction  1,  could 
have  interfered  with  the  rate  coefficient  measurements.  This 
would  have  resulted  in  decreased  kx  values.  However,  since  the 
measured  rate  coefficients  are  independent  of  the  [S02]  ranges 
and  the  initial  [BC1]  concentrations,  cf.  Table  I,  and  since  plots 
of  In  [BC1]  relative  VS  [S02]  are  linear,  we  conclude  that  fluorescence 
from  SO  did  not  significantly  affect  the  measurements.  Fur¬ 
thermore,  BC1  was  monitored  at  four  different  wavelengths,  cf. 
Table  I,  all  of  which  gave  similar  results. 

The  k\  measurements  were  fitted  to  the  expression  k(T)  = 
AT*  zxp{-E/RT)  using  a  three-parameter  nonlinear  regression 
analysis.10  The  value  of  «,  determined  from  the  fit,  was  rounded 
to  one  decimal  place,  and  a  two-parameter  nonlinear  regression 
analysis  was  done.  The  resulting  expression  is 

fc1(470-1690  K)  =  7.7  X  lO'^r/K)30  X 

exp(-1680  K/  T)  cm3  molecule'1  s*1  (3) 
with  variance  and  covariance  matrix  elements11 

{cJA)2  =  0.00142,  a2  =  1 164  K2,  a^A  =  1.193  K 

(4) 

This  method  is  advantageous,  because  it  prevents  round-off  errors 
in  the  value  of  n  in  the  final  fit  and  simplifies  the  error  analysis. 
Combining  the  covariance  matrix  elements  by  the  propagation 
of  errors  technique  yields  ±2 <r*  precision  limits  of  8%  at  470  K, 
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3%  at  1080  K,  and  4%  at  1690  K.  Allowing  10%  uncertainty 
associated  with  the  flow  profile  factor7  and  a  conservative  20% 
for  further  systematic  errors  yields  ±2 ok  confidence  limits  of 
24%  at  470  K  and  23%  at  1080  and  1690  K. 

BC1  +  N20  Reaction.  Preliminary  measurements  showed  a 
dependence  of  k2  on  the  average  flow  velocity,  indicating  that 
continued  BC1  production  in  the  reaction  zone  was  occurring. 
Addition  of  50-500  ppm  of  02  just  downstream  from  the 
microwave  discharge  removed  this  dependence,  and  this  procedure 
was  therefore  adopted.  In  Table  II,  the  rate  coefficient  mea¬ 
surements  on  reaction  2  over  the  temperature  range  670-1000 
K  are  summarized.  Above  about  1000  K  the  measured  rate 
coefficients  were  observed  to  be  dependent  on  pressure,  suggesting 
that  decomposition  of  N20  is  interfering  with  the  measurements. 
Hence  these  are  not  included  in  the  table.  The  k2  measurements 
were  fitted  to  the  Arrhenius  expression  k(T)  =  A  exp (-E/RT) 
to  give 

Jk2(670-1 000  K)  -  5.1  X  10'12  X 

cxp(-6340  K  jT)  cm3  molecule'1  s'1  (5) 
with  matrix  elements 

(aJA)2  *  0.0764,  cE  =  56400  K2,  o^/A  -  65.1  K  (6) 

Combining  these  matrix  elements  by  the  propagation  of  errors 
technique  yields  ±2<r*  precision  limits  of  18%  at  670  K,  7%  at 
830  K,  and  10%  at  1000  K  and  ±2<r*  confidence  limits  of  28% 
at  670  K,  24%  at  830  K,  and  25%  at  1000  K. 

Model-Based  Transition-State  Theory  (MTST) 

Recently,3  we  developed  a  method  for  applying  transition- 
state  theory  to  reactions  of  metal-containing  species.  As  in  CTST , 
the  rate  coefficient  is  given  by12 

k(T)  exp(-A£*  /RT)  (?) 

h  Mb 

where  A£*  is  the  difference  in  zero-point  levels  of  the  transition 
state  and  the  reactants,  q*  is  the  partition  function  for  the  transition 
state,  and  q\  and  q B  are  partition  functions  for  the  reactants.  In 
the  MTST  method,  we  use  molecular  models  to  estimate  the 
vibrational  frequencies  and  rotational  constants  that  are  needed 
to  evaluate  q  * .  For  the  molecular  model,  we  assume  that  harmonic 
restoring  forces  oppose  changes  in  bond  lengths,  angles  between 
valence  bonds,  and  dihedral  angles.  The  potential  energy  V  is 
thus  given  by13 

2V-  £FnAr,2  +  £f9,A e2  +  J/r*Ar*2  (8) 

i  J  k 

where  Fr,  F*,  and  Fr  are  force  constants  for  changes  in  internuclear 
distances  r,  valence  angles  0,  and  dihedral  angles  r  from  their 
respective  equilibrium  values.  As  described  previously,3  the 
MTST  approach  is  suitable  for  reactions  with  transition  states 
which  arise  from  an  early  crossing  of  neutral  and  ionic  potential 
surfaces,  particularly  for  many  oxidation  reactions  of  metal- 
containing  species.  For  such  transition  states,  one  can  make 
reasonably  good  estimates  of  the  geometrical  parameters  and 
force  constants  based  on  Hammond’s  postulate.14  For  example, 
for  an  exothermic  reaction  it  can  be  assumed  that  the  internuclear 
distances,  valence  angles,  and  force  constants  of  the  transition 
states  are  closely  related  to  those  properties  of  the  reactants.  The 
parameters  required  by  the  treatment  that  cannot  be  obtained 
from  properties  of  the  reactants  are  estimated,3  such  that  the 
calculated  vibrational  frequencies  of  transition  states  are  in  good 
agreement  with  ab  initio  computations15”17  on  similar  reactions. 
We  have  shown18  that  for  transition  states  of  eight  reactions,  one 
of  the  reactants  being  a  B-  or  Al-containing  species  in  each  case, 
calculated  vibrational  frequencies  show  an  average  difference 
with  ab  initio  results  of  less  than  about  100  cm-1.  With  such 
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TABLE  I:  Summary  of  Rate  Coefficient  Measurements  for  BC1  +  SO;* 


reation 

zone 

length  P 
(cm)  (mbar) 


cm-3) 


[SO:]* 

(1015 

cm'3) 


v 

(m  r1) 


T 

(K) 


k  ±  fft 

(cm3  molecule-1  s-1) 


reation 

zone 

length 

(cm) 


P 

(mbar) 


[Ml 

(10»7 

cnr3) 


[SO,]* 

(10‘5 

cm'3) 


d 

(m  r1) 


T 

00 


k±  at 

(cm3  molecule"1  rl) 


10 

40  5 

22 

0.6 

25c 

35 

1308 

4.35  ±0.33  (-13)* 

10 

55.3 

3.9 

3.3 

48 

29 

1027 

1.89  ±0.12  (-13) 

VO 

40  6 

2  2 

1  l 

19* 

36 

1316 

4.80  ±0.35  (-13) 

10 

66.3 

4.5 

1.6 

35 

19 

1073 

2^0  ±0.22  (-13) 

10 

64  1 

3  6 

1.2 

55e 

34 

1295 

4.91  ±0.37  (-13) 

20 

46.7 

3.1 

1.1 

34 

27 

1080 

2.00  ±0.18  (-13) 

20 

64  0 

3  6 

0.6 

49* 

34 

1296 

3.93  ±0.32  (-13) 

20 

25.9 

1.6 

0.8 

32 

35 

1138 

2.65  ±0.25  (-13) 

20 

27  9 

1  2 

0.5 

27* 

91 

1622 

1.33  ±0.12  (-12) 

20 

30.5 

2.0 

1.7 

53 

57 

1102 

2.39  ±0.17  (-13) 

10 

27  9 

1  2 

1  0 

38c 

92 

1637 

1.65  ±0.13  (-12) 

20 

39.6 

2.5 

1.1 

32 

35 

1126 

2.46  ±0.26  (-13) 

20 

437 

2  2 

0.4 

21c 

53 

1626 

1.38  ±0.13  (-12) 

10 

39.7 

2.4 

2.6 

56 

37 

1187 

3.61  ±0.27  (-13) 

10 

48  7 

2  2 

0.8 

4CK 

52 

1623 

1.40  ±0.11  (-12) 

10 

89.3 

5.5 

1.4 

45 

21 

1171 

3.02  ±0.18  (-13) 

10 

31  7 

3  2 

5.3 

33f 

18 

723 

3.62  ±0.36  (-14) 

10 

50.3 

3.0 

1.5 

43 

21 

1212 

3.01  ±  0.22  (-13) 

20 

31  6 

3.3 

5.4 

15' 

18 

702 

2.43  ±0.23  (-14) 

10 

32.1 

1.8 

2.3 

66 

42 

1258 

4.49  ±  0.34  (-13) 

20 

82  4 

9  6 

8.0 

27' 

12 

621 

1.13  ±0.08  (-14) 

20 

32.0 

1.8 

1.2 

33 

43 

1277 

3.77  ±  0.30  (-13) 

20 

31  2 

1  4 

1.0 

37c 

101 

1628 

1.11  ±0.09  (-12) 

20 

44.1 

2.6 

1.7 

52 

46 

1226 

3.45  ±0.22  (-13) 

20 

31  6 

1  4 

1.0 

27' 

101 

1651 

1.15  ±0.09  (-12) 

20 

19.9 

1.8 

3.2 

37 

30 

802 

4.99  ±0.51  (-14) 

10 

53  9 

24 

0.8 

24' 

66 

1650 

1.46  ±0.16  (-12) 

—  20 

60.7 

5.8 

2.6 

17 

16 

762 

4.20  ±0.30  (-14) 

10 

62  0 

2  7 

0.7 

8 f 

34 

1685 

1.24  ±0.11  (-12) 

20 

39.7 

3.9 

5.0 

29 

19 

746 

3.33  ±0.23  (-14) 

JO 

29  3 

1  3 

1.3 

17/ 

72 

1691 

1.45  ±0.13  (-12) 

20 

49.0 

5.1 

4.4 

13 

11 

696 

2.32  ±  0.27  (-14) 

29  5 

1  3 

0.7 

2V 

71 

1691 

1.56  ±0.12  (-12) 

10 

49.1 

5.1 

8.7 

15 

11 

692 

2.48  ±  0.16  (-14) 

20 

33  5 

2.7 

3.4 

25T 

28 

903 

8.02  ±0.72  (-14) 

10 

83.2 

9.6 

6.2 

13 

8 

627 

1.41  ±0.11  (-14) 

TO 

67  6 

5.5 

1.1 

17/ 

18 

895 

7.69  ±0.58  (-14) 

20 

28.5 

4.0 

9.6 

11 

10 

523 

3.40  ±0.50  (-15) 

20 

22  5 

1.8 

2.7 

44/ 

55 

904 

1.08  ±0.11  (-13) 

20 

79.8 

12.3 

15.5 

12 

6 

471 

2.12  ±0.19  (-15) 

20 

23  6 

1.8 

2.9 

24 

33 

953 

1.03  ±0.09  (-13) 

20 

63.3 

9.6 

15.9 

19 

6 

476 

2.77  ±0.17  (-15) 

io 

46.1 

3.5 

2.9 

12 

17 

962 

-  68  ±0.78  (-14) 

20 

63.3 

9.1 

15.2 

18 

6 

503 

2.87  ±0.26  (-15) 

10 

28.6 

2.1 

2.6 

21 

37 

977 

1.35  ±0.14  (-13) 

20 

57.7 

8.5 

10.5 

53 

9 

490 

2.31  ±  0.22  (-15) 

10 

44.9 

3.3 

3.2 

20 

30 

975 

1.19  ±0.10  (-13) 

20 

85.4 

12.2 

15.0 

7 

6 

508 

4.05  ±  0.53  (-15) 

10 

78.9 

5.7 

2.8 

14 

17 

999 

1.19  ±0.10  (-13) 

20 

63.9 

8.8 

14.0 

24 

7 

527 

4.65  ±  0.46  (-15) 

10 

27.6 

1.5 

1.5 

14 

66 

1377 

5.99  ±0.91  (-13) 

20 

35.4 

4.7 

1 U 

12 

9 

548 

7.37  ±0.62  (-15) 

10 

85.2 

4.5 

1.0 

16 

29 

1368 

4.62  ±0.48  (-13) 

20 

56.1 

7.2 

12.7 

17 

8 

563 

6.48  ±0.39  (-15) 

10 

54.2 

2.9 

1.2 

24 

40 

1375 

4.91  ±0.49  (-13) 

20 

48.4 

6.1 

10.7 

25 

9 

575 

8.00  ±0.62  (-15) 

10 

54.2 

2.9 

1.2 

4' 

40 

1375 

5.26  ±0.40  (-13) 

20 

54.7 

7.2 

8.4 

38 

11 

552 

4.94  ±  0.60  (-15) 

20 

25.7 

1.4 

0.6 

9* 

82 

1329 

5.23  ±0.69  (-13) 

20 

47.1 

5.7 

8.2 

36 

12 

600 

1.16  ±0.09  (-14) 

20 

25.7 

1.4 

1.2 

39 

82 

1329 

5.53  ±0.54  (-13) 

20 

55.6 

6.9 

7.9 

31 

12 

583 

7.22  ±0.50  (-15) 

20 

59.1 

3.3 

0.5 

76 

42 

1311 

4.89  ±0.39  (-13) 

10 

85.3 

10.4 

11.9 

21 

8 

592 

1.05  ±0.13  (-14) 

20 

33.5 

1.8 

0.6 

88 

50 

1346 

5.37  ±0.45  (-13) 

10 

47.6 

5.2 

7.3 

55 

13 

668 

2.45  ±  0.20  (-14) 

20 

20.1 

1.0 

1.0 

50 

98 

1472 

9.18  ±0.82  (-13) 

10 

77.1 

8.8 

8.7 

34 

11 

635 

1.79  ±0.23  (-14) 

20 

62.2 

3.2 

0.6 

68 

49 

1430 

7.36  ±0.59  (-13) 

10 

40.9 

4.5 

5.6 

93 

17 

657 

2.34  ±0.26  (-14) 

20 

39.5 

1.9 

0.8 

44 

60 

1507 

9.35  ±0.81  (-13) 

20 

20.0 

2.1 

5.3 

84 

18 

674 

1.73  ±0.19  (-14) 

20 

34.5 

1.6 

0.5 

31 

59 

1554 

1.06  ±0.07  (-12) 

20 

24.8 

2.6 

6.1 

21 

16 

678 

1.98  ±0.16  (-14) 

10 

34.6 

1.6 

1.6 

34 

60 

1570 

8.37  ±0.60  (-13) 

20 

20.6 

1.7 

3.4 

37 

28 

897 

8.38  ±0.72  (-14) 

10 

94.6 

4.5 

1.2 

63 

35 

1535 

7.33  ±0.53  (-13) 

20 

51.1 

4.3 

4.3 

51 

22 

854 

5.48  ±  0.33  (-14) 

10 

58.2 

2.6 

0.9 

20 

44 

1602 

1.38  ±0.16  (-12) 

20 

39.5 

3.3 

4.7 

32 

20 

872 

6.12  ±0.48  (-14) 

10 

79.9 

3.8 

0.7 

15 

31 

1536 

1 .09  ±  0.1 1  (-12) 

10 

39.5 

3.3 

4.8 

46 

20 

859 

7.04  ±0.68  (-14) 

10 

19.7 

1.4 

2.5 

31 

38 

1003 

1.69  ±0.19  (-13) 

10 

76.2 

6.7 

3.3 

32 

15 

826 

6.93  ±  0.55  (-14) 

10 

41.1 

2.9 

1.8 

38 

27 

1041 

2.10  ±0.18  (-13) 

a  The  measurements  are  reported  in  the  sequence  in  which  they  were  obtained.  b  In  arbitrary  units 
at  the  272.0-nm  transition. c  BC1  was  monitored  at  the  272.16-nm  transition.  d  Should  be  read  as  (4 
monitored  at  the  266.0-nm  transition,  f  BC1  was  monitored  at  the  272.21-nm  transition. 


Unless  otherwise  indicated,  BC1  was  monitored 
35  ±  0.33)  X  10"13  cm3  molecule1  s~l. €  BC1  was 


estimated  geometrical  parameters  and  force  constants,  we  obtain 
the  vibrational  frequencies  by  a  normal  coordinate  analysis  of 
the  molecular  model,  which  procedure  has  been  described  in  detail 
by  Pauling  and  Wilson.19  Finally,  rotational  constants  are 
calculated  from  the  eigenvalues  of  the  inertia  tensor.1* 

Discussion 

BCI  +  S02  Reaction.  The  thermochemically  accessible 
channels  of  reaction  1 

BCI  +  S02  —  OBC1  +  SO  (la) 

—  B02  +  SC1  (lb) 

are  245  and  26  kJ  mol"1  exothermic,  respectively,  at  0  K.20"22 
Channel  lb  is  probably  inhibited  by  a  substantial  energy  barrier, 
as  all  of  the  reactant  bonds  must  be  broken  and  reformed.  In 
contrast,  we  expect  channel  1  a  to  have  a  slight  barrier  as  only 
the  O-SO  bond  is  broken  in  the  reaction.  In  channel  1  a,  formation 
of  SO (3 2~),  the  ground  state,  is  forbidden  by  the  spin  rule;  however, 
formation  of  SO(J  A)23  is  both  exothermic  (by  169  kJ  mol-1)  and 
spin  allowed.  The  MTST  calculations  which  follow  are  based  on 
the  formation  of  this  state.  Table  III  summarizes  the  rotational 
constants  and  vibrational  frequencies  for  the  transition  state.  The 


first  two  terms  of  eq  7  are  computed  at  several  temperatures  in 
the  range  500-2000  K  and  fitted  to  the  expression  AT.  An 
estimate  of  the  barrier  height  is  made  by  equating  eqs  3  and  7 
at  735  K,  the  midpoint  of  the  Tx  range  investigated,  which  yields 
24.1  kJ  mol"1.  The  present  calculations  lead  to  the  expression 

=  519  X  lO-'V/K)226  X 

exp(-2902  Yk/T)  cm3  molecule*1  s'1  (9) 

Figure  2  compares  the  MTST  calculations  with  the  recommended 
k\(T)  expression,  and  they  are  in  good  agreement  over  the  470- 
1690  K  range.  This  agreement  suggests  that  channel  la,  for 
which  the  MTST  calculations  are  based,  is  the  dominant  pathway. 
In  a  related  paper,3  we  attempted  to  make  rate  coefficient 
measurements  on  the  homologous  A1C1  +  S02  reaction.  We 
failed  to  observe  any  significant  disappearance  of  A1C1  by  reaction 
with  S02,  indicating  that  the  AlCi  reaction  has  a  much  lower 
rate  coefficient  than  the  BCI  +  S02  reaction.  The  lower  rate 
coefficient  is  consistent  with  the  reaction  being  at  least  1 10  kJ 
mol"1  endothermic.20-22 

BO  +  NjO  Reaction.  The  abstraction  pathway  of  reaction  2 

BCI  +  N20  —  OBC1  +  N2  (2a) 

is  631  kJ  mol"1  exothermic  at  0  K.20-22  Similarly  to  our 
investigation  of  the  above  treatment  of  BCI  +  S02,  we  investigate 
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reaction  zone 
length  (cm) 


P  (mbar) 


[M] 

(1017  cur7) 


[N  ;0]  mil 

(1015  cm*3) 


v  (m  s-1) 


k  ±  at 

(cm3  molecule1  r1) 

4.28  ±  0.40  (-16)* 
1.94  ±0.34  (-15) 
1.78  ±0.14  (-15) 
5.50  ±0.73  (-16) 
4.64  ±0.31  (-16) 
6.07  ±  0.59  (-15) 

3.88  ±  0.69  (-15) 

7.14  ±0.64  (-15) 
1.08  ±0.07  (-14) 
1.06  ±0.09  (-14) 

9.14  ±0.97  (-16) 
1.80  ±0.42  (-15) 
1.07  ±0.10  (-15) 
1.48  ±0.12  (-15) 
2.59  ±0.39  (-15) 

2.89  ±0.28  (-15) 
4.01  ±0.43  (-15) 
4.42  ±0.40  (-15) 
7.09  ±  0.60  (-15) 

8.23  ±  0.57  (-15) 
5.56  ±0.43  (-15) 
5.74  ±  0.39  (-15) 
4.29  ±  1.04  (-16) 
5.95  ±  1-06  (—16) 
5.22  ±0.45  (-15) 
4.27  ±  0.52  (-15) 
5.45  ±  0.74  (-15) 


10  4/.1 

•  The  measurements  are  repon^n  the  ^ 

the  272.0- nm  transition. *  Should  be  read  as  (4.28  ±  0.40) 

TABLE  m:  Transition-State  Parameters  for  C1BOSO -  1700 

g^metrical  parameter^  forc^onsunt* - 10-'3  p" 


For  each  measurement,  BC1  was  monitored  at 


MC1B)  1  -72 

KB°)  [9i 

K  OS)  '-43 

KSO)  1  43 

8(ClBO) 

fl(BOS)  l20o 

0(OSO)  120 

r(ClBOS)  0° 

r(BOSO)  90° 

rotational  constants  (B\B2B-i) 
vibrational  frequencies  v( 


force  constant _ _ _ _ 

FAC  IB)  35 

Fr(BO)  -20 

FA  OS)  5.15 

FASO)  10-3 

/XCIBO)  0.2 

F,(BOS)  0.6 

FAOSO)  1-68 

Fr(CIBOS)  0.02 

Fr(BOSO)  0.02 

7.95  X  Kb4  cm"3 
50,  65, 122.  237,  520,  839, 
896, 1300  cm-1 


T  (K) 

1000  700 


T  IQ-1 3 


0  We  assume  trans  arrangements  for  the  atoms  in  the  C1BOS  and 
BOSo'gwups1*  Units  of  r  are  10-10  m.  ‘  Units  of  F, are  102  N  m  units 
of  F>  and  F,  are  Kb18  N  m  rad"2. 

the  kinetics  of  this  reaction  path  using  MTST  with  j^P®^*®***^ 
summarized  in  Table  IV.  The  first  two  terms  of. ar*  “  £ 
in  the  500-2000  K  range  and  fitted  to  the  expression  AT*.  6 
eqS  5  and  7  at  the  mid  T+  point  802  K,  we  back  out  a  barrier 
height  of  53.2  kJ  mob1,  which  combines  with  the  calculated  A 
and  n  values  to  give 
W=7.98Xl0-18(r/K)2-°'x 

exp(-6399  K/T)  cm3  molecule-1  s'1  (10) 

As  shown  in  Figure  2,  eq  10  is  roughly  in  agreement  with  the 
recommended  k2{T)  expression  in  the  670-1000  K  te”Pfra“f* 
range.  We  obtained  similar  results  in  a  previous  study 
homologous  A1C1  +  N20  reaction.  The  two  reactions  are  dose  y 
related  in  that  both  have  fairly  large  activation  barriers  and  both 
reactions  have  exothermic  channels  leading  to  O-atom  abstracts  . 
BC1  +  O2  Reaction.  The  abstraction  reaction 

BCl  +  02  — *■  OBC1  +  O  UD 

is  298  kJ  mob>  exothermic  at  0  K.2^  We  again  ^culate 

transition-state  parameters,  which  are  summanz^inTaWe 

for  the  abstraction  channel.  By  equating  the  MTST  expression 


I  10-'a  F  BCUNjOq 


1000/7  (K-3) 

Figure  2.  Arrbemus  The  Ml 

channels,  as  calculated  by  the  MTST  method. 

to  the  experimental  rate  coefficient  expression1  at  820  K,  we 
calculated  barrier  height  of  31.0  kJ  mob1,  which  yields 

fcMTST=5.17X10-,7(r/K)l'77x 

exp(-3730K/f)  cm3  molecule*  s  (12) 

As  shown  in  Figure  3,  eq  12  is  in  good  agr«ment  wi^  Uie 
experimental  measurements.  In  contrast  to  the  present  study, 
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TABLE  IV:  Transition-State  Parameters  for  C1BONN*  TABLE  VI:  Transition-State  Parameters  for  C1BOCO* 

geomctricai  parameter4^  force  constant  geometncal  parameter4*4 _ force  constant* 

L72  FA CIB)  Ts  riOB)  1.72  FXC1B)  3.5 

HBO)  1.94  Fr(  BO)  -2.0  KBO)  1.94  FABO)  -2.0 

r(ON)  1.18  FA  ON)  5.85  KOC)  1.16  FAOC)  8.0 

r(NN)  1.13  Fr(NN)  18.7  r(CO)  1.16  F,(CO)  16.0 

0(C1BO)  120°  F^CIBO)  0.2  tf(ClBO)  120°  FrfClBO)  0.2 

9(  BON)  90°  F#(BON)  0.6  *(BOC)  90°  FA  BOC)  0.6 

0(ONN)  160°  F#(ONN)  0.67  0(OCO)  160°  F^OCO)  0.79 

r(ClBON)  0°  Fr(ClBON)  0.02  r(ClBOC)  0P  Ft(C1BOC)  0.02 

r(BONN)  0°  Ft(BONN)  0.02  r(BOCO)  0°  Fr(BOCO)  0.02 


rotational  constants  (B&Bi)  2.46  X  10"3  cm"3 

vibrational  frequencies  vt  58,  123,  170,  292,  596,  839, 

996,2221  cm"1 

a  We  assume  trans  arrangements  for  the  atoms  in  the  C1BON  and 
BONN  groups.  *  Units  of  rare  10“10m.  *  Units  of  F,  are  1CPN  nr1;  units 
of  F#  and  Fr  are  10-1*  N  m  rad-2. 


TABLE  V:  Transition-State  Parameters  for  C1BOO* 


geometrical  parameter42' 

force  constant* 

r(ClB) 

1.72 

FA  CIB) 

3.5 

ri  BO) 

1.94 

FABO) 

-2.0 

riOO) 

1.21 

FAOO) 

5.9 

0(C1BO) 

120° 

FKC1BO) 

0.2 

0(BOO) 

120° 

FABOO) 

0.6 

r(ClBOO) 

0° 

Fr(ClBOO) 

0.02 

rotational  constants  (B\BiBs)  1.00  x  10"2  cm-3 
vibrational  frequencies  75,  131,317,841,  1116  cm"1 


a  We  assume  a  trans  arrangements  for  the  atoms  in  the  C1BOO.  b  Units 
of  r  are  10"10  m.  *  Units  of  F,  are  102  N  nr1;  units  of  F#  and  Fr  are  1CH* 
N  m  rad*2. 

T  (K) 

1700  1000  700  500 


1000  '  T  (K-1) 

Figure  3.  Arrhenius  plots  of  measurements  of  the  rate  coefficients  for 
the  BCI  +  O2  (open  squares)  and  BC1  +  CO2  (open  circles)  reactions 
from  previous  studies.1*2  The  full  lines  represent  rate  coefficient 
expressions  for  the  O-atom  abstraction  channels,  as  calculated  by  the 
MTST  method. 

wc  previously  found3  that  the  A1C1  +  02  reaction  MTST 
calculations,  which  were  based  on  an  O-atom  abstraction 
mechanism,  were  inconsistent  with  experiments.  The  A1C1  +  02 
reaction  differs  from  the  analogous  BCI  reaction  in  that  formation 
of  the  abstraction  products  OA1C1  +  O  is  an  endothermic  process 
with  a  barrier3  of  at  least  66  kJ  mol"1.  On  the  basis  of  the  ab 


rotational  constants  (B1B2B3)  2.26  x  10**3  cm"3 

vibrational  frequencies  r(  57, 122,  169,  303,  662,  842, 

1117,2185  cm"1 

a  We  assume  trans  arrangements  for  the  atoms  in  the  C1BOC  and 
BOCO  groups. 4  Units  of  rare  10"l0m. c  Units  of  F,  are  KPN  nr1;  units 
of  F§  and  Fr  are  10"u  N  m  rad"2. 

initio  study  by  Chen,  Hase,  and  Schlegel,17  we  suggested  that  the 
main  channel  of  A1C1  +  O2  reaction  involves  the  formation  of  an 
adduct. 

BO  +  CO2  Reaction.  The  O-atom  abstraction  reaction 

BCI  +  C02  —  OBC1  +  CO  (13) 

is  266  kJ  mol"1  exothermic  at  0  K.20*22  The  calculated  transition- 
state  parameters  for  this  path  are  summarized  in  Table  VI.  By 
equating  the  MTST  expression  to  the  experimental  rate  coefficient 
expression2  at  1080  K,  we  calculate  a  barrier  height  of  73.2  kJ 
mol"1,  which  gives 

W=l-23X10-,7(r/K)i05X 

exp(-$804  K/  T)  cm3  molecule"1  s_1  (14) 

Figure  3  shows  a  comparison  of  eq  14  with  experimental  data 
from  the  previous  study.2  The  calculations  are  in  approximate 
agreement  with  experiment  over  the  770-1830  K  range.  A 
comparison  can  be  made  between  the  BCI  and  A1CI  reactions 
with  C02.  As  is  the  case  for  the  BCI  and  A1C1  reactions  with 
02,  the  BCI  4*  C02  reaction  has  an  exothermic  abstraction  channel 
while  the  channel  leading  to  formation  of  OA1C1  +  CO  has  an 
enthalpy  barrier  of  at  least  98  kJ  mol"1 .  Because  the  abstraction 
channel  for  the  A1C1  +  C02  reaction  is  endothermic,  other 
channels  resulting  in  adduct  formation  can  dominate,  especially 
at  low  temperatures.3 

Conclusions 

Kinetic  measurements  on  BCI  reactions  with  oxygen-containing 
oxidants  have  been  made  over  wide  temperature  ranges.  We 
used  a  TST-based  method  to  calculate  pre-exponentials  for  four 
BCI  reactions  and  obtained  good  agreement  with  experiments. 
We  showed  that  for  these  reactions  rate  coefficients  can  be 
reasonably  well  predicted  over  a  range  of  temperatures,  given  a 
rate  coefficient  measurement  at  a  single  temperature  to  obtain 
the  E  factor.  For  exothermic  processes,  i.e.,  the  N20  reactions, 
A1C1  and  BCI  behave  similarly,  as  may  be  expected  for  cogeners. 
The  observed  differences  for  the  other  reactions  are  in  accord 
with  the  thermochemistry. 
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Activation  barriers  for  series  of  exothermic  homologous  reactions. 

IV.  Comparison  of  measurements  to  theory  for  reactions  of  s>p  atoms 
with  N20 
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The  kinetics  of  the  AI+N,0^A10+N2  reaction  has  been  studied  in  a  high-temperature  fast-flow 
reactor  (HTFFR).  The  "  expression  *(520-1030  K)=4.6X10'12  (T/K)0,5  exp(-778  K/T) 
cm1  molecule'1  s'1  has  been  derived  from  a  nonlinear  regression  analysis  of  the  measured  rate 
coefficients.  2 a  precision  limits  are  about  ±5%  and  accuracy  limits  are  estimated  to  be  about 
±23%.  Combining  the  data  with  a  296  K  literature  measurement  yields  *(296-1030 
K)  =  5.8X  10~15  (T/K)1 37  exp(  —90.4  K/T)  cm3  molecule'1  s'1,  with  2<r estimated  accuracy  limits  of 
about  ±27%.  The  semiempirical  approach  used  previously  to  predict  activation  barriers  for  s'  and 
j2  metal  atom  reactions  with  N,0  has  been  modified  to  allow  predictions  of  group  13  atom  reactions 
with  N-iO.  The  activation  barriers  for  the  B,  Ga,  In.  and  T1  atom  reactions  are  calculated  from  the 
experimental  activation  barrier  of  the  A1  +  N20  reaction.  ©  1995  American  Institute  of  Physics. 


I.  INTRODUCTION 

In  the  first  paper1  of  this  series  it  was  shown  that  the 
experimental  activation  barriers  E ,  and  as  a  result  the  rate 
coefficients  of  reactions  of  ground  state  s'  and  s2  metal  at¬ 
oms  with  N20,  are  closely  related  to  the  sums  of  the  s-p 
promotion  energies  and  ionization  potentials  of  the  atoms.  E 
was  defined  by 

k(T)  =  ATn  exp(  - EIRT)  0) 

which  expression  gave  the  best  correlations  for  n  —  0.5.  A 
semiempirical  treatment  was  developed  to  calculate  E  for  all 
the  s 1  and  s2  metal  atoms  reactions  with  N20,  using  the 
experimentally  determined  activation  barrier  of  the  Cu  +  N20 
reaction.  In  the  second  paper,2  the  treatment  was  refined  for 
5 1  metals  for  which  the  reaction  coordinate  had  to  be  taken 
as  a  variable  instead  of  a  constant.  In  the  third  paper,3  the 
treatment  for  s 2  metals  was  further  modified  and  reactions 
with  02  were  considered.  Here,  the  calculations  are  extended 
to  reactions  of  Z(s2p l)  atoms  with  N20,  where  the  barrier 
heights  for  the  group  are  determined  from  measurements  on 
the  A1  +  N20  reaction,  reported  here. 

II.  EXPERIMENT 
A.  Technique 

The  measurements  on 

A1  +  N,0—  A10  +  N:  (2) 

for  which4  A W298  is  -345  kJ  mol'1,  were  performed  in  an 
HTFFR  (high-temperature  fast-flow  reactor).  This  facility5,6 
and  the  operational  procedures7  have  been  described  in  de- 
tail.  Briefly,  a  vertical  mullite  (McDanel  MV-30)  reaction 
tube  (60  cm  long,  2.2  cm  i  d. )  was  radiatively  heated  to  the 
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desired  temperature  by  columns  of  resistively  heated  SiC 
rods  inside  an  insulated,  water-cooled  vacuum  housing.  A1 
vapor  was  produced  by  resistively  heating  an  Al-wetted  W 
coil  in  a  flow  of  Ar  carrier  gas.  Mixtures  of  N20  in  Ar,  at 
volume  flow  rates  of  between  2%  and  7%  of  the  main  Ar 
flow,  were  injected  into  the  Al/Ar  gas  stream  through  a  mov¬ 
able  inlet.  The  total  molar  flow  rate  of  gas  through  the  inlet 
was  kept  constant  during  a  given  rate  coefficient  measure¬ 
ment.  Further  downstream,  relative  A1  concentrations  were 
determined  by  laser-induced  fluorescence  (LIF)  or  atomic 
resonance  absorption  spectrometry  (ARAS).  ARAS  was  used 
as  a  check  on  the  LIF  results.  These  two  optical  measure¬ 
ment  techniques  have  an  important  difference.  LIF  measures 
the  concentration  of  A1  at  the  axis  of  the  reactor,  whereas 
ARAS  measures  the  average  concentration  over  the  full 
width  of  the  reactor.  Comparison  thus  also  gives  an  indica¬ 
tion  of  the  effect,  if  any,  of  the  flow  profile.  For  the  LIF 
measurements  a  pulsed  Lambda  Physik  EMG  101  excimer/ 
FL2002  dye  laser  in  combination  with  a  KDP  doubling  crys¬ 
tal  was  used.  The  A1  (2S-2P°)  transition  at  265.2  nm  was 
pumped  and  the  fluorescence  intensity  was  measured  through 
a  262  nm  (26  nm  full  width  at  half-maximum)  filter  by  an 
EMI  9813  QA  photomultiplier  tube  (PMT),  connected  to  a 
Data  Precision  Analogic  6000/620  100  MHz  transient  digi¬ 
tizer.  For  the  ARAS  experiments  a  hollow  cathode  lamp  was 
used  to  generate  the  A1  line  radiation  at  309.3  nm  and  the 
transmitted  radiation  intensity  was  monitored  with  a  mono¬ 
chromator  in  combination  with  the  PMT  connected  to  a  Kei- 
thley  614  electrometer. 

Rate  coefficient  measurements  were  made  in  the  station¬ 
ary  inlet  mode7  at  observed  reaction  zone  lengths  of  10  or  20 
cm,  under  pseudo-first-order  conditions  with  A1  as  the  minor 
reactant.  [Al]relalive  is  defined  as  [Al]/[A1], ;  where  [Al],  is  the 
concentration  of  Al  at  the  lowest  [N20]  used.  The  [N20]  was 
typically  varied  by  a  factor  of  5.  In  LIF,  [Al]  relative ^F/  F  (  , 
where  F{  is  the  initial  fluorescence  and  F  the  fluorescence  in 
the  presence  of  N20.  In  ARAS,  the  Lambert- Beer  law  gives 
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[Al]rdal,vc*ln{/n//)*  where  /„  is  the  light  intensity  passing 
through'the  reactor  without  Al  present  and  /  is  the  intensity 
with  Al  atoms.  Rate  coefficients,  k, .  were  obtained  from  the 
slope  of  linear  plots  of  ln[AI]„.,;l„v,.  vs  [N20].  A  weighted 
linear  retzression  was  used  to  determine  kr  and  rr^  lor  each 
measurement.* 

The  gases  used  were  99.998%  Ar  from  the  liquid  Ar 
(Linde),  and  0.53%  N20  (99.99%)  in  Ar  (99.999%)  from 
Matheson,  the  calibration  of  which  was  checked  against  pure 
N20.9 

B.  Results  and  discussion 

The  measured  rate  coefficients  k,  and  the  experimental 
conditions  under  which  they  were  obtained  are  summarized 
in  Table  I.  The  parameters  varied  include:  total  pressure  P, 
corresponding  to  total  concentration  [M ],  initial  fluorescence 
intensity  F  or  initial  percentage  absorption  (both  a  measure 
of  initial  [Al]),  average  gas  velocity  v.  and  reaction  zone 
length.  Our  reported  k,  values  were  determined  to  be  inde¬ 
pendent  of  these  parameters  by  examining  plots  of 
l k(T)-k,\lk(T ).  where  k(T)  is  obtained  from  Eq.  (3)  below 
vs  the  particular  parameter.  For  each  parameter  the  k ,  data 
were  found  to  be  randomly  scattered  about  the  zero  line  of 
the  ordinate.  Figure  I  shows  the  k ;  data  taken.  The  range  of 
temperatures  covered  is  520-1030  K.  The  lower  limit  was 
determined  by  the  heating  of  the  Ar  bath  gas  by  the  Al 
source.  The  upper  limit  was  set  by  the  onset  of  potential 
dissociation  of  N.O.11' 

A  weighted  nonlinear  fit"  of  the  data  to  Eq.  (I)  with 
n=0.5  yields 

k(520-l030  K) 

=  4.6 x  1  0  ”  !  2( 77K)°  5 

Xexp(  -6.469//?r)cm3  molecule'1  s'!,  (3) 

where  E  is  expressed  in  kJ  mol  l.  The  covariance  matrix 
elements1"  are  calculated  to  be  =5.02X10  A  , 
aA£=4.84XlO"3  A£,  and  4  =  4.83 X  I0‘3  E2 .  These  yield 
±2 ak  precision  limits  varying  from  ±2%  to  ±7%.  Fitting 
only  the  LIF  results  would  have  led  to  a  k{T)  expression 
2%-\\%  higher  than  Eq.  (3),  while  considering  only  the 
ARAS  data  yields  an  expression  I0%-13%  lower.  The  later 
set  consists  of  less  data  and  is  therefore  less  precise.  None¬ 
theless.  this  suggests  a  small  systematic  difference,  well 
within  the  ±20%  systematic  errors  and  ±10%  flow  profile 
factor  uncertainty  always  allowed  for  in  this  type  of 
work6,7,1013  Using  these  factors  leads  to  accuracy  intervals 
varying  from  ±22%  to  ±24%.  The  ht  is  shown  in  Fig.  I. 

The  only  previously  reported  kinetic  measurements  on 
the  A!  +  N:0  reaction  resulted  in  a  £<296  K)  of 
<  1 1  ±  1 ) X  10“ 12  cm3  molecule  !s  '.  It  was  obtained  by 
Mitchell  and  co-workers  using  their  static  cell  laser- 
photolysis,  laser-fluorescence  technique.14  Their  uncertainty 
value  pertains  to  precision  only.  Allowing  tor  possible  sys¬ 
tematic  errors  of  similar  magnitude  as  in  our  work,  leads  to 
( 1 1  ±2  to  3) X  10” 12  cm3  molecule' 1  s'1.  Extrapolation  of 
Eq.  (3)  to  296  K  yields  ( 5 .8 ±  1 .8)X  10” 12  cm' 
molecule-1  s-!.  The  joint  uncertainty  limits  thus  do  not  quite 


overlap.  The  value  of  *=0.5  in  Eq.  (3)  was  chosen  for  uni¬ 
formity  with  the  work  on  the  .v1  and  s~  metals,'  "  and  may  not 
be  the  best  choice  here.  Also  extrapolation  over  a  large  range 
of  T  i.e.,  from  ( 1.9  to  3.4) x  10  1  is  a  somewhat  uncertain 
procedure.  There  is  no  reason  to  doubt  either  data  set.  In¬ 
deed,  we  recently13  made  measurements  on  the  Cr  +  N20  re¬ 
action  in  a  MHTP  (metals  high-temperature  photochemistry 
reactor)  from  278  to  1150  K,  which  showed  excellent  agree- 
ment  with  the  static  cell  results  at  298  and  348  K.,16  as  well 
as  with  a  shock  tube  study17  from  1130  to  2570  K.  In  the 
same  work  we  reported  similar  agreement  between  a  MHTP 
study  of  the  Cr+HCI  reaction  from  811  to  1414  K  and  HT- 
FFR  measurements  from  880  to  1449  K.  As  the  various  tech¬ 
niques  for  studies  of  elementary  metal  atom  reactions  are 
showing  good  agreement,  see  also  Refs.  5  and  18,  it  appears 
justified  to  combine  the  two  sets  of  results  for  the  present 
reaction.  Using  the  fitting  procedures  given  above  leads  to 

£(296-1030  K) 

=  5.8X10'I3(77K)1  37 

Xexp( -0.752//?7)cm3  molecule-1  s”1  (4) 

with  4  =  3.19  A",  crAn  =  —0.428  A,  aAE  =  2M  AE, 
4  =  7.86XI0”\  (TnE—  -0.382  £,  and  crE=2.1\  E"  are 
yielding  ±2 crk  precision  limits  from  ±3%  to  ±21%  and  ac¬ 
curacy  limits  from  ±23%  to  ±31%.  Figure  2  shows  that  this 
expression  leads,  in  the  temperature  range  of  the  present 
measurements,  to  a  fit  nearly  as  good  as  that  resulting  from 
Eq.  (3). 


111.  THEORY 

A.  Activation  barriers  of  Z(s2p1)+N20(X  3I+) 
reactions  using  semiempirical  configuration 
interaction  (SECI)  theory 

Following  the  treatment  of  Ref.  1,  the  wave  function  if) 
of  the  activated  complex  (NNOZ)*  can  be  taken  as  the  linear 
superposition  of  three  wave  functions  ,  iA2,  and  corre¬ 
sponding  to  the  properties  of  three  resonating  structures. 
These  result,  respectively,  from  the  covalent  ground-ground 
state  interaction,  a  covalent  interaction  where  a  valence  s 
electron  is  promoted  to  a  p  state,  and  an  ionic  interaction 
between  the  positively  charged  atom  and  negatively  charged 
N20.  Thus, 

?//=  c  |  4/\  +  c i  ifri  +  c ' 

where  ck  (£  =  1.2.3)  are  the  expansion  coefficients.  The 
height  of  the  barrier  is 

£=  J  dr. 

According  to  the  variation  principle,14  the  wave  function 
for  the  system  is  that  for  which  E  is  at  a  minimum.  This 
condition  can  be  expressed  as  set  of  simultaneous  linear 
equations. 

2  ck(Hlk-±fkE)  =  0.  I-  1,2,3.  (7) 

t  - 1 
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TABI.H  1  Summars  <>t  rate  coefficient  measurements  ot  AI-N:0.J 


Reaction 
/one  length 
(cml 

F 

Imbar) 

|,V/ 1 

no17  cm  ') 

[NjOL, 

<10l2cm  ') 

F 

(arbitrary 

units) 

Iniliai  Al 
absorption 

l'/r) 

V 

(ms  ' ) 

7 

<K> 

k  r  <r. 

i cm  molecule  ‘s  ') 

20 

18. 2 

2.1 

3.7 

84 

28 

623 

4.49-0.42  ( - 1  bh 

20 

18.2 

2.1 

3.6 

69 

28 

634 

4.33-0.43  (-11) 

20 

18.2 

2.0 

3.5 

62 

29 

646 

4.33-0.40  (-11) 

20 

34.6 

3.7 

3.8 

52 

27 

672 

4.73 ±0.47  (-11) 

20 

34.6 

3.7 

3.9 

51 

26 

670 

4.35 ±0.40  (-11) 

10 

34.5 

3.9 

4.0 

70 

26 

646 

4.25±0.46  (-11) 

10 

16.6 

1.9 

5.0 

17 

21 

629 

4.75 ±0.65  (-11) 

20 

44.9 

4.4 

2.6 

47 

17 

745 

5  II  ±0.40  (-11) 

20 

44.9  . 

4.2 

2.6 

57 

18 

765 

5.4!  ±0.45  (-11) 

20 

44.7 

4.2 

2.5 

42 

18 

771 

5.91  ±0.46  (-11) 

20 

34.1 

3.9 

2.5 

46 

18 

640 

4.39±0.45  (  —  11) 

20 

34.2 

3.8 

4.1 

72 

19 

655 

4. 1 2 ±0.36  (-11) 

20 

20.3 

12 

2.9 

102 

35 

665 

3.46-0.37  (-11) 

10 

47.9 

4.5 

6.0 

49 

17 

767 

5.29 ±0.39  (-11) 

10 

47.9 

4.5 

6.0 

52 

17 

77! 

5.44±0.45  (-11) 

10 

28.3 

2.6 

4.4 

68 

29 

776 

5.59±0.50  (-11) 

10 

54.4 

5.0 

6.0 

51 

17 

7% 

6.3!  ±0.54  (-11) 

!0 

52.9 

4.4 

4.1 

68 

18 

864 

7.52±0.75  (-11) 

20 

48.7 

3.8 

4.6 

47 

22 

92H 

6.76±0.53  (-1!) 

20 

48.4 

3.8 

3.2 

29 

22 

927 

7. 18±0.48  (-11) 

20 

45.7 

3.6 

3.4 

49 

30 

921 

6.86±0.62  (-11) 

20 

49.3 

3.9 

4.2 

38 

30 

915 

7.28±0.54  (-11) 

20 

19.4 

1.5 

2.6 

51 

40 

938 

6.65±0.75  (-11) 

20 

23.7 

1.8 

3.8 

57 

47 

931 

6.71  ±0.80  (-11) 

20 

22.7 

1.8 

3.0 

40 

60 

921 

6.00±0.62  (-11) 

20 

22.6 

2.1 

3.1 

34 

33 

784 

5.77±0.57  (-11) 

20 

22.6 

2.1 

3.9 

36 

33 

789 

4.22i0.35  (-11) 

20 

25.8 

2.4 

4.5 

31 

34 

792 

5.07 ±0.66  (-11) 

20 

31.0 

2.8 

4.2 

26 

36 

796 

5.01  ±0.51  (-11) 

10 

20.0 

1.9 

4.9 

43 

31 

774 

4.37±0.55  (-11) 

10 

20.0 

1.9 

5.7 

34 

31 

775 

3.72 ±0.44  (-11) 

10 

34.2 

3.2 

7.1 

21 

25 

775 

3.53 ±0.45  (-11) 

20 

21.7 

3.0 

4.7 

37 

27 

516 

2.40±0.26  (-11) 

20 

28.9 

3.5 

5.8 

44 

22 

601 

2.52±0.25  (-11) 

20 

17.3 

1.9 

12.9 

57 

47 

671 

3. 19±0.38  (-11) 

20 

17.3 

1.9 

12.9 

56 

47 

671 

2.93 ±0.32  (-11) 

10 

17.3 

1.9 

13.2 

62 

46 

656 

3.58±0.4i  (-11) 

10 

16.0 

1.8 

13.6 

68 

45 

656 

3.61  ±0.59  (-11) 

10 

20.8 

2.2 

12.2 

63 

50 

675 

4.96 ±0.47  (-11) 

20 

16.1 

1.4 

9.1 

24 

57 

826 

4.57±0.49  (-1!) 

20 

16.1 

1.4 

10.8 

48 

57 

826 

4.56±0.48  (-11) 

20 

14.7 

1.3 

II. 1 

46 

55 

834 

6.05 ±0.66  (-11) 

20 

14.7 

1.3 

10.9 

32 

56 

847 

5.97 ±0.76  (-11) 

10 

14.7 

1.3 

11.4 

34 

54 

813 

4. !7±0.57  (-11) 

10 

14.7 

1.3 

13.1 

39 

54 

812 

4.53 ±0.74  (-11) 

10 

14.7 

1.2 

12.0 

56 

58 

881 

5.95±0.87  (-11) 

10 

14.7 

1.2 

11.8 

47 

60 

896 

6.53  ±0.74  (-11) 

10 

14.8 

1.2 

11.8 

44 

60 

906 

4.94 ±0.66  (-11) 

10 

14.8 

1.2 

11.6 

28 

61 

917 

6.I6±0.78  (-11) 

20 

16.1 

1.3 

9.7 

52 

63 

913 

7.58 ±0.82  (-11) 

20 

16.1 

13 

11.3 

69 

63 

909 

7.22±0.75  (-11) 

20 

19.8 

1.6 

9.3 

35 

66 

901 

7.55 ±0.68  (-11) 

20 

19.9 

1.6 

9.3 

26 

66 

900 

7.94±0.70  (-11) 

20 

26.0 

2.5 

4.6 

9 

26 

743 

3.07±0.24  (-11) 

20 

39.5 

3.8 

4.8 

25 

21 

751 

4.09 ±0.40  (-11) 

20 

39.5 

3.8 

4.8 

22 

21 

752 

3.77 ±0.28  (-11) 

20 

16.4 

1.7 

8.2 

21 

50 

681 

3.91  ±0.41  (-11) 

20 

16.4 

1.7 

10.9 

25 

51 

683 

3.51  ±0.35  (-11) 

10 

16.3 

1.9 

17.5 

24 

47 

632 

3.35 ±0.37  (-11) 

10 

17.2 

n  2 

13.9 

16 

49 

568 

3.01  ±0.34  (-11) 

10 

17.3 

2.2 

27.7 

31 

49 

574 

2.50±0.30  (-11) 

20 

17.2 

2.1 

27.3 

28 

50 

580 

2.31  ±0.24  (-11) 

20 

17.2 

2.1 

16.4 

21 

50 

582 

2.66  ±0.3 1  (-11) 

20 

16.8 

1.3 

8.1 

15 

51 

963 

5.63±0.58  (-11) 

20 

20.3 

1.5 

7.4 

14 

56 

991 

6.36±0.55  (-11) 

20 

20.3 

1.5 

7.4 

15 

56 

992 

6.63 ±0.59  (-11) 

20 

19.2 

1.4 

7.1 

26 

77 

996 

7.7 1  ±0.70  (-11) 

10 

19.5 

1.5 

11.1 

15 

74 

968 

6.03 ±0.65  (-11) 

10 

19.5 

1.5 

18.4 

20 

74 

967 

5.57 ±0.55  (-11) 

20 

32.8 

3.8 

8.4 

16 

19 

627 

3.46 ±0.29  (-11) 

20 

37.9 

2.7 

4.6 

15 

26 

1024 

5.62 ±0.39  (-11) 

20 

37.9 

2.7 

4.6 

17 

26 

1029 

5.75 ±0.47  (-11) 

‘The  measurements  are  reported  in  the  sequence  in  which  they  were  obtained. 
hShould  be  read  as  ( 4.49:!: 0.42 )x  10" 
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T(K) 


10-VTlK'1) 

FIG.  I.  Arrhenius  plot  of  the  rule  coefficient  Jola  of  Al  +  N:0  reaction: 
Q — LIF  used  as  the  monitoring  technique,  - — ARAS  used  as  the  monitor¬ 
ing  technique. 

Here  Hlk  represents  the  interaction-energy  integrals  and  l(k 
the  overlap  integrals,  which  are  defined  as 

Hlk=  j  dr.  A„  =  J  >hf  dr.  (8) 

The  overlap  integrals  are  approximately  equal  to  zero,  A/Jt=0 
for  because  each  ot  the  three  resonating  structures  are 
distinctly  different  from  the  others.  This,  along  with  normal¬ 
ization  of  Eq.  (5),  results  in  \kk  —  J  dr~  1  and 

id= i.  w 

k  =  1 

where  the  ck  represent  the  fraction  of  character  contributed 
by  the  respective  structures  to  the  activated  complex.  Simi¬ 
larly,  Eq.  (7)  reduces  to  the  secular  equation 

//,,-£  H  i:  H]} 

H  ,2  H1:~E  HZ}  =0.  (10) 

Hyy  H2}  Hu~E 

H  ||  is  the  energy  of  the  Zi.s~p 1 1  +  N:0(X  'A  )  interac¬ 
tion.  The  barrier  height  tor  this  structure  cj{s~pl)  equals 
H  |  j  .  For  each  senes  of  homologous  reactions  q  is  constant; 
it  is  the  one  quantity  to  be  determined  from  experiment  in 
these  treatments.  The  second  structure  results  from  the  at¬ 
tractive  covalent  interaction  of  Z is'p")  with  N:0  where  spin 
conservation  requires  N-.0  to  be  in  a  triplet  state,  whereas  tor 
the  5*  and  s~  metals  N-.0  is  in  a  singlet  state  in  the  second 
structure.  The  barrier  to  form  this  structure  includes  the  pro¬ 
motion  energy  PE  of  Z( .v/?1)  to  the  Z(.v'/> "I  state  and  the 
excitation  energy  of  N:0(X  '1  )  to  the  first  N:0(;A  )  state, 
minus  the  amount  by  which  the  interaction  is  attractive.  Fol¬ 
lowing  the  previous  work1-  this  term  is  estimated  to  be 
q{szp 1  )/2.  Thus, 


T(Ki 


lfp/m1) 


FIG.  2.  Arrhenius  plot  of  al!  the  rate  coefficient  data  for  Al  +  N:0  reaction: 
O— Present  work.  □ — Parnis.  Mitchell.  Kamgan.  and  Placket!  t  Ret.  12). 


H22=PE+q(sZp')/2.  OH 

The  third  structure  results  from  the  ionic  interaction  ot 
Z+(.r)  +  N,0"(2n),  which  involves  the  attractive  Coulomb 
and  the  repulsive  Born  potential  between  these  species.  The 
energy  of  the  barrier  resulting  from  this  structure  is  evaluated 
as  before" 

e~  8—  1 

//,,=  IP-EA-  — — r.  <12) 

A0  O 

where  IP  is  the  ionization  potential  of  Z,  EA  is  the  vertical 
electron  affinity  of  N:0,  e  is  the  elementary  charge.  R{)  is  the 
sum  of  the  univalent  radii  of  Z+  and  O  .  and  <5  is  the  average 
Born  exponent  of  Z*  and  O".  Values  for  these  properties  and 
PE  along  with  their  sources  are  given  in  Table  II. 

The  interaction-energy  integrals  are  estimated  to  be 
equal  to  those  for  the  alkali,  alkaline  earth,  and  transition 


TABLE  11.  Physical  properties  used  to  calculate  energy  integrals  tor  group 
13  atom  reactions  with  NT). 


z 

PE 

(kJ  mol  1 1 

IP-EA 
ikJ  mol  [ ) 

Rt} 

HO  "cm) 

S 

B 

34V 

mis1’ 

1  ^  V 

M)J 

Al 

347 

705 

2.53 

7.0 

Ga 

455 

706 

2.5X 

S.O 

in 

4  10 

776 

2.7X 

X.5 

Tl 

54 1 

S07 

3.05 

0.5 

‘The  values  correspond  to  the  lowest  v - p  promotion  energy  obtained  from 
the  atomic  energy  levels  in  Ret.  24. 

"Values  are  obtained  from  the  ionization  potentials  in  Ret.  25  and  the  ver¬ 
bal  electron  affinity  of  NT).  -217  kJ  mol  which  corresponds  to  the 
NT)  negative-ion-siatc  II  from  Ret  26. 

Univalent  radii  are  derived  from  the  effective  ionic  radii  in  Ret.  27  The 
univalent  radius  ot  O  is  taken  to  be  I  .76X  10  cm. 

■•‘The  average  Born  exponent  equals  (df+7)/2  where  the  value  ot  7  is  the 
Born  exponent  ot  O  and  the  Si  value  is  that  ot  Z  taken  trom  Ret.  -8. 
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TABLE  II!  Comparison  of  experimental  and  calculated  activation  barriers  for  group  13  atom  reactions  with 
NLO 


z 

(fraction) 

C  2 

(fraction) 

c] 

(fraction) 

A 

(cm1  molecule-1  s'1] 

)  (kJ  mol  1 

^e»pi 

)  (kJ  mol " 1 ! 

Tjngc 

i  (K) 

B 

0.85 

0.09 

0.06 

8.4X10-12 

21.3 

22.  lb 

300 

Ref.  20J 

A) 

0.79 

0.09 

0.12 

4.6X  10" 12j 

6.5 

6.5 

520-1030 

This  work 

Ga 

0.8!  ‘ 

0.06 

0.13 

5.9X10’12 

14.3 

I4.4e 

296 

Ref.  21J 

In 

0.81 

0.07 

0.12 

6.3XKT12 

13.4 

T! 

0.87 

0.04 

0.09 

7.1  x  10"12 

26.7 

■■Experimental  value  of  the  A)  +  N20  reaction:  A  coefficients  for  all  other  reactions  are  calculated  using  simple 
collision  theory. 

hTh\s  £„ p,  value  is  obtained  from  the  300  K  rate  coefficient  and  the  A  coefficient  of  column  4. 

"This  £cxpl  value  is  obtained  from  the  296  K  rate  coefficient  and  the  A  coefficient  of  column  4. 

Reference  from  which  the  experimental  data  were  obtained. 


metal  atom  reactions.2  Thus  H  !2  and  Hn  are  equal  to  fi  and 

is  equal  to  fill ,  where  fi=~  107.9  kJ  mol  1 . 

B.  Results  and  discussion 

The  value  of  q(s2p[)  is  chosen  such  that  the  height  of 
the  barrier  £  of  Eq.  (10)  equals  the  experimental  activation 
barrier  £expt  of  the  A1  +  N20  reaction  (Sec.  II  B).  From  Eq. 
(3)  this  yields  q(s2p 1 )  =85.7  kJ  mol'1.  Using  this  value,  and 
Eqs.  (7)  and  (9),  the  activation  barriers  £  and  the  fraction  of 
character  contributed  by  each  structure  can  now  be  calcu¬ 
lated  for  the  other  boron  group  atoms.  The  results  of  these 
calculations  are  given  in  Table  III.  Table  III  also  shows  the 
comparison  of  the  calculated  activation  barriers  to  the  de¬ 
rived  experimental  activation  barriers  £expt.  For  B  and  Ga, 
these  experimental  activation  barriers  are  determined  by  set¬ 
ting  the  rate  coefficients  found  in  room  temperature 
studies20,21  equal  to  AT0,5  exp (-E/RT)  in  which  the  A  are 


IP  +  PE  (kJ  mol1) 


FIG  3.  Plot  of  the  experimental  and  calculated  activation  barriers  i crosses) 
of  Z *N:0— ZO-*-N:  reactions  vs  the  sums  of  the  \~p  promotion  energies 
and  ionization  potentials  of  Z.  Alkali  metal  (open  triangles),  alkaline  earth, 
and  transition  metal  (open  circles),  and  boron  group  (open  squares)  atom 
reactions  are  distinguished. 


taken  as  the  values  given  by  simple  collision  theory22 
Rl(%irkB/fi)0'5.  Here,  kB  is  the  Boltzmann  constant  and  p, 
the  reduced  mass.  The  predicted  bamers  for  B  and  Ga  are 
within  0.7  kJ  mol-1  of  the  derived  experimental  values,  us¬ 
ing  the  lowest  s-p  promotion  energies.  Using  instead  the 
activation  barrier  from  the  combined  experimental  rate  coef¬ 
ficient  expression,  Eq.  (4),  yields  values  within  1.2  kJ  mol-1 
of  the  predicted  values.  This  confirms  that  varying  the  value 
of  n  has  little  influence  on  the  accuracy  of  the  predictions. 
The  £exp(  for  In  and  T1+N20  could  not  be  estimated  because 
no  experimental  data  are  available  for  these  reactions. 

In  the  present  treatment.  Sec.  Ill  A,  we  assumed  the 
lowest  s-p  promotion.  The  activation  barriers  were  also 
computed  using  the  lowest  p-s  promotion  energies  and  the 
average  of  the  promotion  energies  [(j-p )  +  {/?-.? )]/2.  In  all 
cases  the  calculated  activation  barriers  were  too  high  for  bo¬ 
ron  and  too  low  for  gallium,  indicating  that  the  lowest  p-s 
promotion  energies  do  not  correlate  with  the  activation  bar¬ 
rier  of  these  atoms. 

In  Fig.  3  the  experimental  and  calculated  activation  bar¬ 
riers  of  the  present  work  are  compared  to  those  for  the  s 1  and 
s 2  metals  atom  N:0  reactions.1 2  As  can  be  seen,  good  agree¬ 
ment  between  theory  and  experiment  is  obtained  for  all  three 
groups.23  Each  group  shows  a  linear  increase  of  the  activa¬ 
tion  barrier  with  the  sums  of  the  ionization  potentials  and 
s-p  promotion  energies,  which  properties  determine  the  en¬ 
ergies  of  the  resonating  structures,  cf.  Eqs.  (11)  and  ( 12).  As 
the  magnitude  of  IP  +  PE  increases,  the  resonance  interaction 
between  these  structures  and  the  activated  complex  weakens, 
thereby  increasing  the  barrier  £. 

In  summary,  the  semiempirical  treatment,  applied  to  r1 
and  s2  metal  atom  reactions  with  N20,  has  been  extended  to 
s2p{  atoms  by  modifications  to  the  treatment  and  use  of  the 
experimental  barrier  from  one  s2p]  atom  reaction,  namely 
Al  +  N-,0.  Comparison  to  the  experimental  data  available  in 
the  literature  indicates  good  agreement  between  theory  and 
experiment. 
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APPENDIX  F 

The  AIO  +  02  Reaction  System  Over  A  Wide  Temperature  Range. 

David  P.  Belyung  and  Arthur  Fontijn* 

High-Temperature  Reaction  Kinetics  Laboratory,  The  Hermann  Department  of  Chemical 
Engineering.  Rensselaer  Polytechnic  Institute,  Troy,  NY  12180-3590 

(Received 

The  kinetics  of  the  AIO  +  02  -action  system  has  been  studied  in  a  high-temperature  fast-flow 
reactor  (HTFFR)  from  7  to  100  mbar.  The  relative  AIO  concentrations  were  monitored  by  laser- 
induced  fluorescence  at  the  B2£-X*X  and  CTl-X^X  transitions.  Below  about  1000  K  the  AIO 
consumption  is  pressure-dependent,  i.e.  leading  initially  to  AlOj.  with  log  k(305-1010  K)  = 

-25.36  -  1.69  log  (T  /  K)  cm*  molecule-2  sr>,  as  obtained  by  weighted  linear  regression.  At  higher 
temperatures  the  system  becomes  pressure-independent  with  AIO,  and  O  as  the  likely  products.  For 
this  process  k(  1200-1690  K)  =  7.7  x  1  O' 10  exp  (-10008  K/T)cm3  molecule-1  S'1.  The  accuracy 
estimates  are  discussed  in  the  test.  The  reactions  are  suggested  to  proceed  through  an  OAIOO  addition 
complex  of  similar  structure  as  identified  in  solid  matrices. 


INTRODUCTION 


Reactions  of  aluminum  species  play  a  role  in  many  high-temperature  environments.1"4  We 
have  studied  a  large  number  of  A 1,  A1C1  and  AlO  reactions  over  wide  temperature  ranges,  by  use  of 
the  HTFFR  (high-temperature  fast-flow  reactor)  technique.5"8  Two  early  studies  of  the 
AlO  +  O2  system,  at  respectively  1400  and  300  K,9-10  suggested  a  temperature-independent  rate 
coefficient,  attributed  to  abstraction 

AlO  +  O2  — ^  AIO2  +0  ( 1 ) 

This  interpretation,  which  indicates  no  activation  energy,  is  in  conflict  with  some  subsequent 
information,  such  as  the  currently  recommended  endothermicity 1 1  for  reaction  ( 1),  A H°  =  96  ±  40 
kJ  mol-1.  Recently,  we  showed  that  several  reactions,  which  were  considered  in  the  literature  to 
proceed  by  simple  abstraction,  i.e.  AlO  +  CH4  and  Cr  +  O2,  to  be  more  complex.6-12  The  earlier 
AlO  +  O2  measurements  could  have  similarly  hidden  an  addition  process,  as  they  showed 
considerable  scatter.  Since  then  the  HTFFR  technique  has  been  much  improved.5-12  A  re¬ 
investigation  of  the  AIO/O2  system  thus  was  indicated.  We  demonstrate  here  that  there  are  in  fact 
two  reactions,  reaction  (1)  which  under  the  present  conditions  dominates  above  1200  K,  and  an 
addition  reaction 


AlO  +  O2  +  Ar  — >  AIO3  +  Ar 
which  dominates  below  1000  K. 


(2) 


TECHNIQUE 

The  HTFFR  facility  used  in  this  work  and  the  data  handling  procedures  have  been  previously 
described.5-14  Briefly,  a  vertical  mullite  (McDanel  MV-30)  reaction  tube  (60  cm  long,  2.2  cm  i.d.), 
was  heated  to  the  desired  temperature  by  resistively  heated  SiC  rods  inside  an  insulated,  water-cooled 


vacuum  housing.  AlO  was  produced  by  several  methods.  At  room  temperature  and  for  a  few  checks  at 
547  K,  production  was  achieved  by  passing  Ar  with  approximately  100  ppm  of  trimethyl  aluminum 
(TMA)  and  1  ppm  of  02  through  a  microwave  discharge.  This  technique  was  limited  to  reaction  zone 
pressures  from  7  to  28  mbar  and  could  not  be  used  at  higher  temperatures.  For  the  other  experiments 
A1  vapor  was  produced  by  resistively  heating  an  Al-wetted  W  coil  in  a  flow  of  Ar  carrier  gas.  A  trace  of 
oxidizer,  0.5%  N,0  or  5%  02  in  Ar,  at  a  4  x  10-3%  0f  the  main  Ar  flow,  was  passed  through  a  Pt 

tube  to  a  location  just  downstream  from  the  vapor  source  to  rapidly  convert  A1  to  AlO5'81516  This 
production  technique  limited  observations  at  547  K  to  P  <  32  mtar  due  to  die  condensation  of  A1  and 
AlO  on  the  reactor  walls,  but  at  higher  temperatures  allowed  observations  up  to  100  mbar,  i.e.  over  the 
full  normal  HTFFR  pressure  range.  For  the  T  £  850  K  experiments  only  the  5%  O j  was  used. 

Downstream  from  the  AlO  production  zone,  0,/Ar  mixtures  at  flow  rates  between  1%  and  8% 
of  the  main  Ar  flow  rate  were  introduced  through  a  movable  inlet  at  10  or  20  cm  upstream  from  the 
observatton  windows.  There,  relative  AlO  concentrations  were  measured  by  laser-induced  fluorescence 
(LIF).  A  pulsed  Lambda  Physik  EMG  101  excimer  /  FL  2002  dye  laser  in  combination  with  a  KDP 
doubling  crystal  was  employed.  Just  before  the  reactor  window,  the  laser  beam  passed  through  two 
light  baffles.  13  made  of  metallic  concentric  brass  cones  with  4  and  2.5  mm  diam.  exits  respectively,  to 
reduce  scattered  light  in  the  reactor.  Two  transitions,  B2I-X2I  and  C2n-X2I  were  used.  For  the 
former.  AlO  is  pumped  on  the  464.8  nm  (1,0)  band  and  the  fluorescence  is  observed  through  a  482nm 
(20  nm  full  width  at  half  maximum)  interference  filter,  i.e.  mainly  at  the  486.6  nm  (1,1)  band.  The 
B^-X’-I  system  could  only  be  used  up  to  1200  K  because  of  interference  by  background  rad.ation 
from  the  reactor  walls.  For  measurements  above  1200  K  and  for  a  comparison  near  826  K  the  weaker 

LIF  from  the  C2n-X2I  (0.0)  transition  at  302.2  nm 1 8  was  used  in  combination  with  a  30 1  nm  ( 1 1  nm 

full  width  at  half  maximum)  interference  filter.  The  fluorescence  was  detected  by  an  EMI  98 1 3Q  A 
photomultiplier  tube,  connected  to  a  Data  Precision  Analogic  6000/620  100  MHz  transient  digitizer. 
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The  gases  used  were  Ar  (99.998%)  from  the  liquid  (Linde),  O,  (99.6%)  from  Linde,  0; 

(99.98%)  from  Matheson  and  5%  02  (99.99%)  in  Ar  (99.995%)  from  Scott.  The  rate  coefficients 
obtained  with  these  O2  sources  were  indistinguishable.  The  5%  O2  and  0.5%  02  (99.99%)  in  Ar 
(99.995%)  from  Matheson,  0.5%  N20  (99.99%)  in  Ar  (99.995%)  from  Matheson,  and  TMA  (98%) 
from  the  liquid  TMA  (Johnson  Matthey)  were  used  for  AlO  production. 

Measurements  were  made  under  pseudo-first-order  conditions,  [AlO]  « [02],  in  the  stationary 
inlet  mode.19  Rate  coefficients  kj,  in  second  order  units,  for  a  given  temperature,  pressure,  and  average 
velocity  were  obtained  by  using  five  02  concentrations  providing  variation  by  a  factor  of  about  5.  kj 
and  <7  iq  were  calculated  by  applying  a  weighted  linear  regression20  to  plots  of  In  [A10]relative  vs.  [O2]. 
These  yielded  straight  lines  with  slopes  equal  to  -kjt,  where  t  is  reaction  time. 

RESULTS 

The  measurements  are  divided  into  two  sets.  In  one,  the  measurements  were  made  at  varying 
pressures  at  temperatures  near  305,  547, 826.  or  1010  K  to  evaluate  the  temperature  dependence  of  the 
pressure-dependent  reaction.  In  the  other,  T  >  1200  K,  where  pressure-independent  rate  coeffictents 
were  obtained,  the  temperature  was  randomly  varied  up  to  1690  K. 

305-1010  K  MEASUREMENTS 

The  kj  measurements  obtained  in  this  range,  and  the  conditions  under  which  they  were  obtained, 
are  assembled  in  Table  I  and  show  a  pressure  dependence.  At  each  temperature,  to  determine  if  the 
data  were  dependent  on  parameters  other  than  [M],  plots  of  the  residuals  [kj  -  k(M)]/k(M)  versus  an 
individual  parameter  were  examined.  Here,  k(M)  is  the  fitting  expression  for  that  particular  temperature, 
eqs.  3,  4,  or  5  below.  These  plots  showed  that  the  third-order  rate  coefficients  are  independent  of  initial 
fluorescence  intensity  F  (a  measure  of  initial  [AlO] ),  average  gas  velocity  v,  and  observed  reaction  zone 

length. 
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There  were  slight  variations  in  temperature  for  each  group  of  measurements.  The  standard 
deviation  of  these  were  computed*'  for  each  data  set  and  combmed,  by  a  propagation  of  etror 
method,22  with  an  estimated  maximum  aT  =  ±  2%  T  for  systematic  errors.  The  resulting  standard 

deviations  for  each  average  temperature  are  shown  in  Table  II. 

The  305  K  measurements  are  near  the  lugh-pressure  limit  as  evidenced  by  the  levelmg  off  of  the 
k,  vs.  [M]  plot.  Fig.  1.  These  data  were  fitted  using  Marquardt's  method,22  with  a  fall-off  expression 

from  T  roe:-4 


log  k(M)  =  log 


k2.0[M] 


l  +  k2,0[M]/k2. 


log  Fc 


l  +  [log(k2i0[M]/k2,oo)f 


(3) 


Here  k2  „  is  the  thtrd-order  rate  coefficten,  for  reaction  (2),  k2._  is  the  high-pressure  second-order  rate 

coefficient,  and  Fc  is  a  broadening  factor.  Eq.  (3)  y.elds  Km  -  5.1  *  '°'30™6  molecule-2  s  ', 
k  =  1. 1  X  10-'2  cm3  molecule-'  s-'.  and  Fc  =  0.752.  To  estimate  Ok2,0  a  nonlinear  fit2  ot  the  data 

was  made  with  only  the  first  term  of  eq.  3.  From  this  the  elements  of  the  covariance  matrix  were 

calculated-5  and<7k2  0  was  found  to  be  0.34  k2.0,  Table  II. 

At  547  and  826  K  linear  k,  vs.  [M]  plots  were  obtained.  Figs  2  and  3.  The  ki  were  fitted  with 
an  expression  purely  proportional  to  pressure. 


k(M)=  k2.o  [M] 


(4) 


From  the  fits22  the  covariance  matrices  were  calculated22  and  the  dk2.„  obtained.  Table  II.  The  547  K 
measurements  allowed  comparison  of  the  various  AlO  production  techniques.  As  shown  in  Fig  2  no 
influence  of  the  method  used  is  evident.  At  826  K  higher  pressures  could  be  studied  than  at  547  K.  due 
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to  the  improved  efficiency  of  the  production  technique.  The  larger  [AlO]  at  this  temperature  also  made  it 
convenient  to  compare  the  LIF  results  from  the  B-—  transition  to  those  of  the  C -I"!  X-X 
transition.  This  is  done  in  Fig.  3,  which  shows  the  independence  of  the  transition  used. 

The  data  at  1010  K  display  linear  behavior  with  an  intercept  indicative  of  a  pressure-independent 
channel.  Fig.  4.  These  data  were  fitted  with  a  simple  linear  expression 

k(M)  =  ki  +  k2,o  [M]  (5) 

using  a  linear  least-squares  algorithm.20  In  this  eq.  k\  is  the  second-order  rate  coefficient.  This 
algorithm  also  yielded  0^2,0  °ki- 

The  temperature  dependence  of  the  third-order  rate  coefficients  is  shown  in  Fig.  5.  The  rate 
coefficients  decrease  with  increasing  temperature,  typical  behavior  for  an  association  reaction.  The  data 
were  fitted  with  a  log  k2,o  =  log  A2,o  +  n  log(T  /K)  expression  using  a  linear  least  squares  fit20  with  the 
transformation  crT  =  2.303  Tojog  T  which  yields 


log  k2.o  (305  -  1010  K)  =  -25.36  -  1.69  log  (T  /  K)  cm6  molecule'2  S'1.  (6) 

The  variance  and  covariance  matrix  elements25  are  cr2i0gA2.o  =  1-45,  oiOgA2.0  n  =  -0.508.  and 
<j-n  =  0. 178.  Using  these  variances  and  covariance  the  2ok20  precision  limits  are  calculated  to  vary 

from  ±76%  at  305  K,  decreasing  to  a  minimum  of  ±22%  at  700  K,  and  increasing  again  to  ±37%  at 
1010  K.  Using  a  ±10%  uncertainty  in  the  reactor  flow  profile,19-26  along  with  a  conservative  ±20% 
for  possible  systematic  errors,  yields  confidence  intervals  varying  from  ±79%  to  ±31%. 
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1200  -  1690  K  MEASUREMENTS 

The  1200  K  to  1690  K  data  are  assembled  in  Table  m.  Residual  plots  showed  uo  dependence  on 
pressure,  initial  fluorescence  intensity  F,  average  gas  velocity  v.  and  mac, ion  rone  length.  As  a 
precaution,  in  determining  the  pammeters  of  the  k,(T)  expresston  below,  only  measurements  at  low 
enough  pressures  to  insure  less  than  20%  contribution  from  the  extrapolaied  tfurd-order  reaction  eq.  6 
were  used.  These  measurements  are  shown  in  Fig.  6  and  show  a  strong  positive  temperature  dependence 
with  no  significant  curvature.  Applying  a  linear  regression  fit? 23 to 
ki(T)  =  Aexp(  -E  K / T)  yields 


k i ( 1 200  -  1690  K)  =  7.7  x  10'10  exp  (-10008  K/  T)  cm3  molecule  1  s  l. 


(7) 


A  nearly  identical  expression  was  obtained  by  (i)  fitting  all  the  data  and  (ii)  by  subtracting  the 
termolecular  component  from  all  the  data  and  fining  the  results.  In  both  cases  the  activation  energy  was 
within  6%  of  the  value  shown  in  eq.  7.  Equatron  7  has  vartances  and  covariance*  of  <7,2  =  0.182  AA 
aAE  =  2.54  x  102  A,  a„d  a£2  =  3.59  x  lO*.  which  yields  precision  limits  varying  from  ±9%  to  ±17% 
depending  upon  temperature.  The  result, ng  estimated  accuracy  limit  is  ±26%  when  the  same  systematic 
etrors  as  above  are  considered.  The  k,  from  the  1010  K  measurements  (above)  was  not  used  as  the 
associated  uncertainly  was  considered  to  be  too  large. 


DISCUSSION 

The  results  of  die  negative  temperature  dependence  of  reaction  (2)  and  the  positive  dependence 


of  reaction  ( 1)  is  that  the  AlO  consumption  rate 


coefficients,  calculated  on  a  second  order  basis,  pass 


through  a  minimum  which  occurs  a,  around  800  -  1000  K  a,  the  pressures  used.  Comparison  of  the  k, 
data  a,  300  K  (Table  I)  to  those  at  1400  K  (Table  III)  shows  little  difference.  Thus  the  temperature 


independence  indicated  for  reaction  ( 1 )  in  the  previous 


work10  was  based  on  a  coincidence.  The  actual 


measurements  from  the  present  and  previous' 


9.10  work  are  in  reasonable  agreement  as  can  be  seen  in 
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Figs.  5  and  6.  For  the  comparison  in  Fig.  5  the  individual  second-order  rate  coefficients  from  ref.  10 
were  divided  by  the  pressure  and  averaged.  There  now  is  also  no  conflict  between  the  endothermicity 
for  reaction  (1)  of  96  ±  40  kJ  mol-1  and  the  activation  energy  of  83  ±  10  kJ  mol-1  obtained  from  eq.  7. 


The  type  of  results  obtained  here  can  be  fitted  to  a  general  scheme  of  the  kind: 


AIO  +  O2 

1 

1 


*■ - *  AI02  +  0 


AI03 


(8) 


,  where  the  dotted  lines  indicate  the  less  likely  channels.  The  path  leading  to  AIO3  appears  a  straight 
forward  termolecular  mechanism.  The  chemically  activated  intermediate  A103  is  also  likely  to  be  the 

precursor  of  the  AIO2  +  O  reaction  products  of  the  high-temperature  second-order  process,  as  the  pre¬ 
exponential  of  the  ki(T)  expression,  7.7  x  lO'10  cm3  molecule-*  s-l,  appears  rather  too  large  for  a 
diatomic  radical  -  diatomic  radical  direct  abstraction  reaction.27  Certainty  with  regard  to  the  mechanism 
of  the  relative  importance  of  direct  and  indirect  bimolecular  processes  is,  however,  generally  difficult  to 

establish.28 

The  large  pre-exponential  of  kj(T)  also  suggests  that  the  AIO3  complex  once  formed  will  not 
significantly  re-dissociate  to  the  reactants  AIO  and  O2.  This  in  contrast  to  the  exothermic  iso-electronic 
BO  +  O2  reaction  system  for  which  no  pressure  dependence  and  a  negative  activation  energy  was 
found.29  This  suggests29-30  a  reaction  involving  an  activated  complex  which  is  not  significantly 
pressure  stabilized.  The  mechanism  for  that  reaction  can  be  summarized  by 

BO  +  O2  BO3  BO2  +  O  ^ 


There  have  apparently  not  been  any  direct  gas-phase  observations  of  A103  species.  However, 

the  A1  +  O2  and  A1  +  O3  reacttons  when  earned  out  in  solid  matrices  have  been  shown  to  yield  a 

product^  *’32  OAIO2  *a.  is  also  the  likely  product  for  AlO,  Oj  addition.  The  mamx  work  further 
suggested  that  a  major  isomer  of  this  stnrcture  is  a  branched  ring  compound  O  -Al^  This  ,s  the  same 

configuration  as  an  SCF/CI  calculation  has  indicated  for  the  B03  intermediate  of  eq.  9.-9  AlO  has  been 
observed  to  form  gaseous  addition  complexes  with  other  molecules.  Thus,  chemiluminescense 
observations  in  the  visible  on  the  A1  +  C02  +  03  reaction  system  have  given  evidence  for  the  formauon 
of  OAICO2 .33  Also,  room  temperature  kinedc  measurement  have  found  drat  to  reactions  of  AlO  with 
C02  and  C2H4  have  pressure-dependent  pseudo-second-order  rate  coefficients.  *  These  C02  and 
C2H4  observations  show  an  interesting  difference  with  the  present  AlO  +  O2  findings,  in  that  toy 
display  a  non-zero  intercept,  indicative  of  a  contribution  of  second-order  process  at  room  temperature. 
That  work  was  done  at  bath  gas  N20  pressures  of  270  to  800  mbar.  Since  in  the  present  work  at 
pressures  an  order  of  magnitude  lower,  with  a  less  efficient  bath  gas,  Ar,  the  second-order  component 
did  not  become  evident  till  about  1000  K  this  may  suggest  a  rather  stable  OAIO2  adduct. 
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TABLE  I 

Summary  of  Rale  Coefficient  Measurements  on  the  AiO  +  02  reaction  taken  below  1050  K 
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TABLE  II 


Summary  of  Third-Oixier  Rate  Coefficients  for  the  AlO/02/Ar  Reaction  System. 


T  ±  Oj 

(K) 

Na 

P  range 
(mbar) 

[M]  range 
(1017  cm'3) 

k2,0  ±  Ok2,0 
(cm^  molecule'2  s'O 

305  ±6 

18 

6.8  -  27.9 

1.6-  6.6 

5.1  ±  1.7(-30)b 

547  ±  19 

14 

7.2  -  32.4 

0.9  -  4.5 

9.1  ±  0.5(-31) 

826  ±24 

15 

11.1  -  102.6 

1.0-  8.9 

5.7  ±  0.3(-31) 

1010  ±26 

30 

11.3  -  81.7 

0.8  -  5.8 

3.0  ±  0.4(-31) 

a)  N  is  the  number  of  ki  measurements  at  this  temperature. 

b)  Should  be  read  as  (5.1  ±  1.7)  x  lO30  cm6  molecule'2  s'1 


TABLE  III  .  10nnkr 

Summary  of  Rate  Coefficient  Measurements  on  the  AlO  +  02  reaction  taken  above  12(X 
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FIGURE  CAPTIONS 

FIG.  1.  Plot  of  kj  vs.  [M]  for  the  AlO/CVAr  reaction  system  at  T  =  305  K 
( O )  individual  rate  coefficients 
( - )  Troe  fit  of  measurements,  eq.  3 

FIG.  2.  Plot  of  kj  vs.  [M]  for  the  AlO/CVAr  reaction  system  at  T  =  547  K 
( O )  vaporization  of  A1  and  5%  O2  used  to  produce  AlO 
( + )  vaporization  of  A1  and  0.5%  N2O  used  to  produce  AlO 
( □ )  microwave  discharge  of  TMA  and  0.5%  O2  used  to  produce  AlO 
(  - )  Proportional  fit  of  measurements,  eq.  4 

FIG.  3.  Plot  of  kj  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  826  K 
( O  )  Measurements  taken  using  the  C2n-X2I  transition  with  LEF 
( □  )  Measurements  taken  using  the  B2£-X2£  transition  with  LEF 
(  -  )  Proportional  fit  of  measurements,  eq.  4 

FIG.  4.  Plot  of  kj  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  1010  K 
( O)  Measurements  taken  using  O2  (99.98%)  from  Matheson 
( □  )  Measurements  taken  using  O2  (99.6%)  from  Linde 
(  -  )  Linear  fit  of  measurements  eq.  5 

FIG.  5.  Plot  of  the  third-order  rate  coefficient  data  of  the  AlO  +  Ot  reaction.  Uncertainties  shown  are 
±2cr  precision. 


(A)  Current  study 
(  -  )  Best  fit  to  current  study,  eq.  6 
(•)  Calculated  from  ref.  10 


FIG.  6.  Arrhenius  plot  of  the  second-order  rate  coefficient  data  of  the  AlO  +  02  reaction. 
(A)  Individual  rate  coefficients 
(-)  Arrhenius  fit  to  data,  eq.  7 
(•)  Recommendation  from  ref.  9 
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